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ABSTRACT 

We provide a theoretical description of synchrotron fluctuations arising from magnetic turbulence. 
We derive an expression that relates the correlation of synchrotron fluctuations for an arbitrary in- 
dex of relativistic electrons to the correlations arising from a particular 7 = 2 index that provides 
synchrotron emissivity proportional to the squared intensity of perpendicular to the line of sight com- 
ponent of magnetic field. We provide a detailed study of the statistics in the latter case assuming 
that the underlying magnetic turbulence is axisymmetric. We obtain general relations valid for an 
arbitrary model of magnetic axisymmetric turbulence and analyze the relations for the particular 
example of magnetic turbulence that is supported by numerical simulations. We predict that the 
synchrotron intensity fluctuations are anisotropic with larger correlation present along the direction 
of magnetic field. This anisotropy is dominated by the quadrupole component with the ratio between 
quadrupole and monopolc parts being sensitive to the compressibility of underlying turbulence. Our 
work opens avenues for quantitative studies of magnetic turbulence in our galaxy and beyond using 
synchrotron emission. It also outlines the directions of how synchrotron foreground emission can be 
separated from cosmological signal, i.e. from CMB or highly redshifted HI emission. For the sake of 
completeness we also provide the expressions for the synchrotron polarization (Stocks parameters and 
their combinations) for the model of axisymmetric magnetic turbulence. 
Subject headings: turbulence ~ ISM: general, structure - MHD ~ radio continuum: ISM. 



1. INTRODUCTION 

Synchrotron fiuctuations both carry information about magnetic fields and interfere with the measurements of cosmic 
microwave emission (CMB) as well as with attempts to measure enigmatic emission from atomic hydrogen (HI) in the 
early Universe. Therefore it is essential to have proper description of the statistics of synchrotron fluctuations. This 
paper concentrates on the statistical properties of the synchrotron intensities arising from turbulent magnetic flelds. 

Galactic and extragalactic synchrotron emission arises from the interaction of astrophysical magnetic flelds and 
cosmic rays (see |Ginzburg 1981 ). In terms of CMB and high redshift HI studies, the most important are galactic 
synchrotron emission, emission provides the largest range of scales for studyin g magnetic fields. Magnetic fields 
are turbulent as observations testify that turbulence is ubiquitous in astrophysics (Armstrong et al. 1995 Lazarian 



2009[ Chepurnov & Lazarian 2010). As relativistic electrons are present in most cases, the turbulence results in 
synchrotron fluctuations, which may provide detailed information about magnetic flelds at different scales, but, at the 
same time, i mpede measures of C MB and high redshift HI. The latter has become a dir ection of intensive discussion 



recently (see |Loeb fc Zaldarr iaga" 2004} [Pen et al.||2008| |Loeb fc Wyithe|l2008[ |Liu et al.|[2009[ ). 

Present and tuture big telescopes, eTg^ Square Kilometer Array (SKA), LOw I'requency ARray (LOFAR), are 
designed to map synchrotron emission with unprecedented resolution and sensitivity. Therefore it is absolutely crucial 
that the information available with these instruments is being used to get better insight fundamental astrophysical 
processes. 

Properties of turbulent magnetic fields that infiucnces most of the processes in diffuse astrophysical environments 
are very much desired through the corresponding studies. The obstacle for such studies is that the relation between 
the underlying turbulence and the synchrotron fiuctuations is not trivial because the synchrotron intensity depends 
both on magnetic field H and on the power index a of relativistic electrons, which emit radiation spiraling in magnetic 
field. The physics of emission is rather non-trivial which results in the expression of synchrotron intensity which are 
proportional to magnetic field in a fractional power, iJ'*', 7 = 1/2(q; + 1) that has been an impediment for a correct 
description of synchrotron fiuctuations. 

For instance, studies of synchrotron intensity fiuctuations in the Milky Way (see Cho & Lazarian 2010 and references 
therein) reveal power law spectra that can be related to magnetic turbulence. How are these spectra related to the 
underlying spectra of magnetic fiuctuations? What is the effect of the electron spectral index on the result? Are the 
fluctuations of magnetic field related to compressions or have Alfvenic nature? Are there additional information apart 
from the magnetic spectrum one can obtain from synchrotron fluctuations? These are the questions that we are going 
to address in the paper. 
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A number of attempts to describe synchrotron fluctuations has been att empted assuming that 7 = 2, which provides 
an approximation to the actual index of galactic synchrotron emission (see Getmantsev|1959 Chibisov fc Ptuskin|1981 



Lazarian fc Shutenkov 
T559; Chepurnov||1998 



1990 



Lazarian fc Chibisov 1991 Chepurnov 1998[ ). Some of these papers (e.g. |Getmantsev 
Che purnov| |1998p did not consider mean magnetic held while others dealt with the spectrum of fluctuations 
that arise from unrealistic idealized isotropic magnetic turbulence. As we discuss further, such an isotropic model 
contradicts the current understanding of magnetic turbulence. 

An attempt to take anisotropy into account was made in Lazarian fc Shutenkov (1990), where the anisotropy of the 
synchrotron fluctuations was analyzed and the quadruple component of the spatial decomposition of the synchrotron 
correlations was used to study the spectrum of synchrotron fluctuations. The work used the explicit integral equations 
solving the inverse problem for the magnetic fluctuations, but the anisotropy was introduced in the problem by 
assuming that the magnetic field can be presented as a superposition of a regular and isotropic random magnetic 
fields. The latter assumption was not supported by the later numerical and theoretical research. 

In short, the limitations of all earlier studies were related both to the incorrect models of astrophysical turbulence 
and to the constrains arising from doing calculation just for a single cosmic ray electron index 7 = 2. The magnetic 
turbulence is known to be anisotropic and the actual index may substantially differs from that single chosen value and 
it is not clear how accurate the description of synchrotron fluctuations stays for 7 not equal to 2. 

The goal of this paper is to address both of the deficiencies of the earlier studies. First of all, recent years have 
been marked by important advance s in the theory of compressible MHD turbulence a nd its statistical description 
e.gjG oldreich & Sridhar 1995; Laz arian fc Vishniac||l9 99 ; Cho fc Vishniac 2000; Maron fc Goldreich|[200T| [Lithwick 

Kowal fc Lazariani|201l)|. I'his understanding of 



fc 



'oldreichi|2001[ iCho et al.l|2002i |Cho fc Lazarian||2002( |2003[ il 
turbulence should be included in tne description of synchrotron fluctuations. In addition, it is absolutely vital to 
understand what is happening with the synchrotron statistics for 7 different from the single value of 7 = 2 that was 
the focus of the earlier studies. 

Our paper relays synchrotron fluctuations with realistic magnetic turbulence for arbitrary electron power law index 
extending the reliability of magnetic held studies. We also propose new ways of studying magnetic flclds, in particular, 
separating contributions from compressible and incompressible magnetic fluctuations. 

In this paper we do not directly address the synchrotron polarization studies. We, however, should mention that the 
existing studies attempting to answer important questions about magnetic fields, e.g. question of magnetic helicity, 
using synchrotron (see Waelkens, Schekochihin fc Enfilin 2009, Junklewitz fc Enfilin 2011) suffer the from same limi- 
tations as the earlier statistical studies of synchrotron intensity, i.e. they also assume 7 = 2 and isotropy of magnetic 
turbulence. We believe that our approach should be useful for these studies. To illustrate our point we address the 
problem of Faraday rotation measure fluctuations arising from anisotropic turbulence. 

Our present paper is complementary to our work on probing turbulence using spectral measurements. In terms 
of those studie s focused on recovering the spectrum of velocit y turbulence from Doppler broadened absorption and 



emission lines (Lazarian fc Pogosyan 2000 2004 2006 



2008 ) the present paper is intended to provide additional 



information about magnetic field fluctuations, 'i'his study is intended to provide the basis for the future development 
of techniques that recover turbulent statistics from observations using not only intensity, but other Stocks parameters 
of synchrotron emission. 

Our study is intended to provide foundations for studying magnetic turbulence in the halo of Milky Way, in supernova 
remnants, in external galaxies as well as in clusters of galaxies (see EnBlin et al. 2010, 2009, Enfilin fc Vogt 2006). As 
the resolution of telescopes gets higher detailed studies of turbulence in astrophysical jets, synchrotron emitting lobes 
get possible. 

In §2 we justify the model of magnetic turbulence that we adopt, in §3 we describe the fluctuations of synchrotron 
emissivity and obtain general expressions for synchrotron fluctuations for an arbitrary spectral index of relativistic 
electrons. In §4 we analyse the anisotropy of synchrotron fluctuations arising from an arbitrary model of axisymmetric 
turbulence, while in §5 we show how this general treatment gets modifled when a model of MHD turbulence that follows 
from numerical simulations is used. We discuss our results in §6 and provide our summary in §7. Appendixes are 
important part of these work. Apart from containing derivations related to the synchrotron intensity correlations they 
also contain the statistics of correlations of other Stocks parameters and useful combinations of the Stocks parameters 
(see Appendix E). 

2. STATISTICAL DESCRIPTION OF MHD TURBULENCE: PERSPECTIVE OF AN OBSERVER 



MH D Turbulence plays a crucial role for the proc esses of cosmic ray propa gation (see 
2011 |, star formation (see [Elmcgrccn fc Scalo 2004] |McKee fc Ostriker||2007|), heat trans' 



Schlickeiser 2002 



er m magnetizec 



Longair 
plasmas 



(see|Narayan fc Medvedev 2001 j |Lazarian||2006| ), magnetic reconnection (see Lazarian fc Vishniac 1999 Kowal et al. 
2009p . Statistical description allows to reveal regular features within chaotic picture of turbulent fluctuations. 'I'he 
famous Kolmogorov description of incompressible hydrodynamic turbulence provides a vivid example of how turbulence 
complexity can be reduced to a simple formula E{k) ~ k^^^, where E{k) is an energy spectrum of turbulent motions. 
MHD tur bulence is more complex than the hydrodynamical one. Magn etic fleld defines the chosen direction of 



anisotropy (Montgomery fc Turner 1981 Shebalin et al. 1983 



Higdonjl984| . For small scale motions this is true even 
in the absence of the mean magnetic field in the system. In this situation the magnetic field of large eddies defines the 
direction of anisotropy for smaller eddies. This observation brings us to the notion of local system of reference, which 
is one of the major pillars of th e modern theory of MHD turbulence. Therefore a correct formulation of the theory 
requires wavelet description (see Kowal fc Lazariaii||2010 ). Indeed, a customary description of anisotropic turbulence 
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using parallel and perpendicular wavenumbers assumes that the direction is fixed in space. In this situation, however, 
the turbulence loses its universality in th e sense that, for instance, t he critical balance condition of the widely accepted 
model incompressible MHD turbulence ( Goldreich fc Sridhar||1995 henceforth GS95) which expresses the equality of 
the time of wave transfer along the magnetic held lines and tne eddy turnover time is not satisfied"'^. 

MHD turbulence is a developing field with its ongoing debates^. However, we feel that among all the existing 
models the pS95 provides the best correspond ence to the existing numerical and observational data (see Beresnyak| 



& Lazarian 2010 Chepurnov & Lazarian 2010). A limitation of the original GS95 model is that it assumes that the 
injection oi energy happens with Vl equal to the Alfven velocity Va- In astrophysical situations turbulence can be 
injected with both Vl < Va and Vl > Va- For the first case the turbulence starts initially in the regime of weak 
turbulence when Alfvenic perturbations weakly interact with each other. As scales of motions decrease, the strength 
of interactions, paradoxically, increases and the turbulence gets stronger. This regime was described analytically in 
Lazarian & Vishniac (1999, henceforth LV99). The regime Vl > Va is the regime of super-Alfvenic turbulence, when 
turbulence initially follows the Kolmogorov cascade and magnetic fiel ds are not dynamically important. At smaller 
scales the turbulence gets into the strong MHD turbulence regime (see Lazarian] [2006 ) . We should add that the exact 
form of the magnetic correlatio n tensor does not foll ow from the GS95 model and should be obtained from numerical 
experiments. This was done in Cho et al. (2002) and Cho & Lazarian ( 2003) . 

The local system of reference is not accessible to an observer who deals with projection of magnetic fields from 
the volume to the pictorial plane. The projection effects inevitably mask the actual direction of magnetic field within 
individual eddies along the line of sight. As a result, the scale dependent anisotropy predicted in the GS95, which taking 
into our considerations above (see footnote 1) should be written as Ay ~ L^^^X^/^, where Ay and A_l are the parallel 
and perpendicular scale of the eddy, respectively, is not valid for the observer measuring parallel and perpendicular 
scales of projected and averaged along the line of sight eddies. As the observer maps the projected magnetic field in 
the global reference, e.g. system of reference of the mean field, the anisotropy of eddies becomes scale-independent and 
the degree of anisotropy gets determined by the anisotropy of the largest eddies which projections are mapped. This 



property of eddies being scale indepen dent in the global system of reference is discussed in the earlier works (e.g. Cho 
et al. ||2002[ [Esquivel fc Lazarian||2005| 



I'he testm g of how ideas of incompressib le MHD turbulence are applica ble in the presence of compressibility was 



performed in |Cho fc Lazarian ( 2002[ 20031 (also see Cho fc Lazarian||2005 for a review). There, a decomposition of 
motions into basic MHD modes (Alfvenic, slow and last) was perforrned to show th at Alfvenic and slow modes keep 
the scaling of the incompressible MHD, while fast modes exhibit isotropy (see GS95, Lithwick fc Goldreich||200T| . 

In a number of situations magnetic turbulence becomes isotropic for the observer. One of them is when the turbulence 
is super-Alfvenic, i.e. the velocity at the injection scale Vl is larger than the Alfven speed Va- Even in the opposite 
case, Va > Vl, the observer may still detect isotropic turbulence if the line of sig ht is directed along t he mean field in 
the volu me. Super-Alfvenic turbulence may take place for molecular clouds (see Padoan et al. 2004) and clusters of 
galaxies ( Brunetti fc Lazarian||2007(@| ) . 

In the presence of the mean magnetic field in the volume under study, an observer will see anisotropic turbulence, 
where statistical properties of magnetic field differ in the directions orthogonal and parallel to the mean ma gnetic field 



The description of a xisymmetric turbulence was given by Batchelor ( 1946 ), Chandrasekhar (11950 ) and later Matthaeus 



fc Smith (1981) and Oughton et al. (1997). In our case it is natural to identify the axis of symmetry with the mean 



(1) 



(2) 



magnetic held and then the index-symmetric part of the correlation tensor can be presented in the following form: 

where the separation vector r = xi — X2 has the magnitude r and the direction specified by the unit vector f . The 
direction of the symmetry axis set by the mean magnetic field is given by the unit vector A and fj, = f • A. The magnetic 
field correlation tensor may also have antisymmetric, helical part, however it does not contribute to synchrotron 
correlations. Hence, further on we shall consider only the index-symmetric part of {Hi{xx)Hj{x2)) - 
The structure function of the field has the same representation 

i((iJ,(xi) - i7,;(x2)) (i?j(xi) - i/j(x2))) = A{r,fi)hfj + B{r,n)S^j + C{r,n)X^Xj + D{r,fi) (^hX^ 

with coefficients A(r, ji) = Aj(0, ^) — A(^{r,^) , . . . etc. In case of power-law spectra, one can use either the correlation 
function or the structure function, depending on the spectral slope. In this case the structure function coefficients can 
be thought of as renormalized correlation coefficients. 

^ The original GS95 model was formulated in terms of global syste m of r e ference and us ed closur e relations which are only valid in this 
syste m of reference. This point was corrected in further publications ( ^Lazarian &: Vishniac| 1999 Ch o fc Vishniac|2000|[Maron fc Goldreich| 
|2001| |Cho et al |2002| l where the necessity of using local reference system detmed by magnetic held on the scale ot turbulence study was 
suggested, justitied arid tested. When in the literature the critical balance and the relation between parallel and perpendicular scales are 
still expressed in terms of parallel and perpendicular wayenumbers, the corresponding wavenumbers should be understood as shorthand 
notations of the inverse scales calculated in the local reference frame, rather than wavenumbers in the system of reference related to the 
mean magnetic field. 

^ Recent debates, for instance, were centered on the role of dynamical alignment or polarization intermittency that could change the 
slope of MHD turbulence from the Kolmogorov slope predicted in the GS95 to a more shallow slope observed in numerical simulations (see 
[Boldy rev 2005, 2006, Bercsnyak fc Lazarian 2006). More recent studies (see |Bercsnyak fc Lazarian 2009 201(5| indicate that numerical 
simulafions may not have enough dynamical range to test the actual spectrum ot turbulence and the llattening of the spectrum measured 
in the numerical simulations is expected due to MHD turbulence being less local than its hydro counterpart. 
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In Appendix |B] we discuss the relation of the correlation tensor and spectrum tensor representation of the axisym- 
metric turbulent fields. An observer studying synchrotron does not map projected magnetic fluctuations directly, as 
the synchrotron emission is a non-linear function of the magnetic field component perpendicular to the line of sight. 
This is the subject of the next section. 

3. 3D SPATIAL CORRELATIONS OF SYNCHROTRON EMISSIVITY 

The goal of this section is to show that one can introduce measures of synchrotron correlation which very weakly 
depend on the electron spectral distribution, i.e. on 7. Whenever possible we provide the prove analytically for selected 
cases. We also show results of numerical calculations. 



The emission has been discussed 
Careful study of the formation 



3.1. Synchrotron correlations 

Synchrotron emission arises from relativist ic electron spiralin g about magnetic fields, 
in many monographs (see Pa cholczyk 1970, Fleishman |2008 and references therein) 
of the synchrotron signal (see Westfold||1959 1 revealed that the signal is essentially non-linear in the magnetic field H 
(B). The origin of nonlinearity is m relativistic effects. Nonlinearity comes from the fact that the signal is formed 
only over the narrow fraction of the electron cycle, and the two leading orders in the deviation of the trajectory 
from the straight line give the same contribution in terms o f I/7. Summat ion of the result over "flashes" is also 
non-straightforward, and produces spread frequency spectrum ( Westfold|1959 1. Situation is remarkably different from 
cyclotron (non relativistic) emission where emission is monochromatic and has intensity just quadratic in the magnetic 
field. 

If the distribution of relativistic electrons is 



then the observer sees intensity of the synchrotron emission is, thus, 



sync 



(X) oc J dzHl{-K) 



(3) 



(4) 



where X — {x, y) is the 2D position vector on the sky and H± — H'i + H"^ is the magnitude of the magnetic field 

perpendicular to the line-of-sight z. Note that 7 = ^(q; + 1) is, generally, a fractional power. 
The relativistic electron power law index a changes from object to object and also varies with energ y of electrons 



For galactic radio halo tested at meter wavelengths, observations indicate that a sa 2.7 (see Pohl 



and therefore 7 ma y vary due to processes of acceleration and loses. In the simplest models 01 shoe 



see 



|Long; 



1996|. The index a 
t acceleration a = 2 



air 



2011[ ) and in the acceleration in turbulent reconnection a = 2.5 (see de Gouveia dal Pino & Lazarian 
2005[ ). I'hose should produce 7 = 1.5 and 7 = 1.75, respectively. In reality, the acceleration of particles is a more 
sophisticat ed process which in case of a shock includes th e formation of the precursor and its interaction with the 



media (see Beresnyak et al.|2009 Malkov &: Diamond||2009| as well as various feedback processes. As a result of these 
and propagation ettects 7 will vary. In this paper we consider variations of 7 in the range from 1 to 4, which covers 
the astrophysically important cases we are aware of. 
In observations, one measures the correlation function of the synchrotron intensity 



Csj/nc(R') — {Isy7ic0^l)Isync{^2)) OC 



X2)) 



(5) 



by averaging over an ensemble of the pairs of the sky measurements at fixed two dimensional separation R = Xi — X2. 
This function is the projection of the three-dimensional correlation of emissivity 



^^.(R,z)-(i/2(xi)i?I(x2)) 

which for homogeneous turbulence depends only on xi — X2 = (R, z). 
The structure function is formally related to the correlation one 



DsyncCR.) = (^(IsyncO^l) " Isync0^2)f'^ OC 2 J dzi j dz (^^^7(0,2) - ^H^i^, z] 



(6) 



(7) 



but it can often be defined even when the correlation function itself has a divergent behaviour. 

Thus, we are presented with the problem of describing the correlation between the fractional powers of the magnitudes 
of the orthogonal projections of a vector field, S,h1 ■ Our main aim is to demonstrate that one can develop statistical 
measures that are insensitive to the parameter 7. Using these statistics simplifies the interpretation of the observational 
signal in terms of the properties of the magnetized turbulence. 

Our approach is the following. Let us note that however complex the statistics of the magnetic field is, the observed 
correlations are functions of just three functional quantities, £,xx{'^) = {Hx{xi)Hx{x2)) , £,yy{^) = {Hy{xx)Hy{x2)) and 
^j.j^(r) = {Hx{'x.-i)Hy(x2)) ■ To study the 7 dependence of the correlation of emissivity ^//t, we shall consider it in the 

space spanned by ^xx,S.yy,S.xy as independent variables, {Cxx,^yy,(.xy)- 



3.2. Statistically isotropic magnetic field and synchrotron correlations 

First we consider the case when statistics of the magnetic field is isotropic, which may, for instance, correspond to 
the super- Alfvenic turbulence, i.e. for the turbulence with the injection velocity much in excess of the Alfvenic one. 
The structure tensor of a Gaussian isotropic vector field, a special case of Eq. ([2]), is usually written in the form 

((F,(xi) - i7,(x2)) (i/j(xi) - i/,(x2))) = {Dll - Dr,N) hfj + D^n^^j , (8) 

where DLL{r) and £)jvAf(^) are structure functions that describe, respectively, the correlation of the vector components 
parallel and orthogonal to point separation r. In case of solenoidal vector field, in particular the magnetic field, two 
structure functions arc related by 

^Dll = --{Dll-Dnn) (9) 
dr r 

which in the regime of the power-law behaviour D^n on r'" leads to both functions being proportional to each other 

To describe our system we chose the x-axis to coincide with the direction of the projected separation R. This can 
be done for analysis of isotropic statistics since there is no other preferred direction. In this case, we find the following 
two correlations 

di = (^[H,{Xi, zi)-H,{X2,Z2)f) = {Dll -Dnn) sin'' + Dnn (10) 
[Hy(Xi,zi)~Hy{X2,Z2)f)^DNN (11) 
[i/,(Xi,zi)-i7,(X2,Z2)]')=0 (12) 

to determine completely the correlations of the orthogonal components of the magnetic field. Here 9 is the angle 
between the 3D separation vector r and the line of sight. 
For convenience we shall also introduce the normalized correlation coefficients 

ci = (i/,(Xi,zi)i7,(X2,Z2))/(if.(Xi,zi)2) = 1- ^di (13) 

C2 = {Hy{Xi,zi)Hy{X2,Z2))/{Hy{Xi,zi)^) - 1 - ^^2 (14) 

which vary from unity for coincident points to zero at large separation. Normalized structure functions, conversely, 
are zero at small separations and reach asymptotically the value of two where the correlation vanishes. 
Let us now investigate the normalized correlations 



Chi = [{Hi (xi) Hi (X2)) - {Hi (x))^ j / [{^Hl i^r) {Hi (x))^ j (15) 

DHi^2(l-iHl) (16) 
where the mean and the variance can be trivially expressed through Gamma functions 

{Hl{x))=2^/'r[l + ^/2]H^ (17) 

'Hl{xf)^2^r[l+j]H'^ (18) 



where H is the amplitude of the magnetic field. Both the mean {Hi (x)) and the variance (^Hl (x)"^^ — {Hi (x))^ 

that determine the average intensity and the rms fluctuations of the synchrotron emission are strong functions of 7 
index. However, as we shall demonstrate, the normalized correlation is not, thus providing a direct measure to the 
correlation properties of the magnetic field. 

Below we identify the normalized correlations of synchrotron emission as the measures which depend on the index 
7 only weakly. The goal is to reduce the problem of relating magnetic field and synchrotron statistics for an arbitrary 
7 to the mathematically trackable case of 7 = 2. 

3.2.1. Special cases of spectral index: "f — 2 and 4 

In case of even integer 7, among which of the most interest are cases 7 = 2 and 7 = 4, the correlation of the 
emissivity can be found analytically to yield 

= I {ci' + C2^) 7 = 2 (19) 



ini = I {ci' + C2') + ^ (ci^ + C2^) + ^^c,W 7 = 4 (20) 

Starting at unity when ci = C2 = 1, the correlations show remarkably little difference throughout the whole range of 
ci and C2, although ^ l'?//^ as Ci — ?> 0, C2 — > 0. 
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In terms of the structure functions 

Dh-i ^ {di + d2) 



{dl 



dl) 



DH-^ = l{di 



rfa) [dl 

' 20 ^ 



d2Y 



7 = 2 
7 = 4 



(21) 
(22) 



Although a different index leads to a change of the amplitude at small separation where linear term dominates by a 
factor 6/5 it does not affect the slope, and as we shall see numerically the difference in the transition to non-linear 
terms will be difficult to observe. 

3.2.2. Case of equal correlation coefficients ci = C2 

Another special test case that can be treated analytically is the case of equal correlation between normal components 
and correlation of longitudinal components. In this case, for arbitrary 7 



i-r[i + 7]r [1 + 7/2]"' 



(23) 



where 2F1 [—2] is the hypergeometric function (see Abramovitz fc Stegun|l965 1. As the Figure [T] shows, this result is 
also nearly 7 independent for 7 e (1.5, 3), and with only minor differences showing only for 7 > 3 





Fig. 1. 



0.2 0.4 

Correlation coefficients , 



0.6 0.8 1.0^^ 0.2 0.4 0.6 0.8 1.0 

'•h'' ("^i) coincident ci = C2 for several 7 £ (1.5,2.8) (left) and 7 £ (2.5,4) (right). 



Ci=C2 



3.2.3. Conjecture 0/7 independence 

Based on the test cases above we conjecture that the normalized correlation function of the synchrotron is insensitive 
to the 7-parameter in the range of our interest. To test this conjecture we compute numerically ^ji/t- for non-integer 
7's. The results presented in Figure [2] confirm insensitivity to 7 in all the range of interest. The maximal deviation in 




Fig. 2. — Left: Example comparison of the correlation coefficients £,^1 for 7 = 2 (top surface) and 7 = 3.5 (bottom surface). Right: The 

difference between two surfaces in the left panel. The systematic separation between those two quite distinct cases is below 5%. Random 
features at 1-2% level are due to limitations of numerical integration. 



correlations are limited to 3% for 7 in the wide range 7 S (1.2,3), while staying within 10% for 7 = 4. 
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Consequently, we conclude that the isotropic correlation of 3D synchrotron emissivity can be well approximated by 
a simple quadratic result of 7 = 2 

IhI^\{c^+c^) (24) 

Dhi - (di +d2)~\ [dl + dl) (25) 

In other words, while the intensity of emission changes appreciably, the correlation properties of the normalized 
intensities depend on 7 only weakly. It is important to note that while the mapping from the magnetic field to 
synchrotron emissivity statistics is essentially non-linear in terms of correlation function, it is linear at small scales 
if the structure functions are used. That makes a structure function the preferable statistics to apply. The effects 
of 7 manifests itself in an amplitude shift of the small scale scaling, an effect difficult to observer without very wide 
dynamical range. 

3.3. Example of statistically anisotropic signal 

While correlation properties of synchrotron emission in statistically isotropic magnetic field are described with two 
functions, in general, full description requires three function D^x^ Dyy and D^y. Study of emission from axisymmetric 
turbulence belongs to this general case. 

Here we shall consider the synchrotron emission in the limit when the anisotropy can be expected to play the largest 
role. This is the limit when a) the symmetry axis is orthogonal to the line-of-sight and b) the field has no variations in 
the direction of the symmetry axis. If we chose x-direction to coincide with the axis of symmetry (z-direction is along 
the line of sight) in this limit 

(^xx = 0, cTxy — 0, ayy = a 

Dxx^O, Dxy^O, Dyy = a^d, (26) 

i.e. the correlations are described just by a single function d. Such a case has physical significance when the axis of 
symmetry is associated with a regular component of magnetic field, so realistic model should include non-zero mean 
field in x-direction, (Hx) ^ 0. Perpendicular to the line of sight component of the magnetic field now has a random, 
y and regular, x components, — Hy A- {Hx)'^- 
The 3D emissivity of the synchrotron is correlated according to 

= {{Hliv,) + {Hx?y" (H^yir,) + {Hx)Y') ({H^y + {Hx?)') (27) 
which is readily computed for 7 = 2 and 7 = 4 

IhI=c\ 7 = 2 (28) 
2 f 3(f74 - c2) \ 



where c is a normalized correlation function of the magnetic field fluctuations perpendicular to the regular component 
of magnetic field. We find that the correlation properties of the synchrotron do not depend at all on the amplitude of 
the regular magnetic component if 7 = 2. The left panel in Figure [3] demonstrates for 7 = 4 that the difference from 
simple 7 = 2 result decreases quickly with the strength of the regular field. Between these two cases, the correlation 
differ by no m ore that 2% along all the range of values when {Hx) > 3(t. Similarly, for all intermediate values of 7, 



7 = 2 formula ( 28 ) provide the limiting behaviour when regular magnetic component is increased. 

Thus, we conclude that 7 dependence of normalized correlations is very weak for anisotropic turbulence in the 
presence of the notable regular magnetic field component. With good precision the result can be approximated by 
7 = 2 formula = {Hy{v^)Hy{v^))\ 

The strongest sensitivity to 7 is expected when the regular component is small, {Hx) ^ cr. While the case of 
strictly one-dimensional {Hx) — 0, {Hx)"^ = turbulent magnetic field is probably of little astrophysical importance, 
it provides a useful, mathematically tractable limit, of strongly anisotropic turbulence. In this limit 

^^^V^rlM, ^,„-^.^^^^,^:iMA ,30) 

= ^^^^ 2^ 2^2^"J ^ (31) 

- V^r[i+7]r[i±^]^^^i 

V^r[i + 7]-rm^ 

The right panel in Figure [3] shows how mapping from the magnetic field to emissivity correlations behave in this limit. 
The 7 = 2 amplitude-independent curve provides the upper bound for both 7 < 2 and 7 > 2 models. 
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Fig. 3. — Left: correlation coefficients (c) for 7 = 4 model in strongly anisotropic turbulence with varying strength of the magnitude 
of the regular component of the magnetic field (dashed curves starting for the lowest one at A = 0). The correlation approaches the 
7 = 2 (solid) curve as (Hx) increases beyond the rms of the random component (Hy)^^'^. Right: 7-dependence in the limit of anisotropic 
turbulence with absent mean field A = 0. The upper curve is 7 = 2 case, same as on the right, the dashed curve corresponds to 7 = 1 and 
the bottom solid one to 7 = 4. 



3.4. 2D correlations of synchrotron intensity in isotropic turbulence and 7 dependence 

Integration of the synchrotron emissivity along the hne of sight provides the observer with the intensity of the 
emission. To demonstrate that it does not modify our conjecture of weak sensitivity to 7, we perform a detailed study 
of 7 = 2 and 7 = 4 isotropic models. 

Line-of-sight projection inevitably includes correlation between widely separated regions of space. Thus one needs a 
statistical model of the magnetic field that is self-consistent at large separations. In physical environment the scaling 
of the turbulence regime can not extend to arbitrary large scales and is expected to saturate at the scales associated 
with energy injection r/ or the size of the turbulent region. We shall consider a simple prescription for the structure 
function of normal components 

DatjvW = D(oo) ( -^^) (33) 



r-j 



which follows the power-law scaling cx r™ at small scales and saturates at large separations at the value equal twice 
the variance of the field component D{oo) = 2 (Hf). To satisfy the solinoidality condition of Eq. ^ the longitudinal 
structure function is 

„ , , D(oo) { r \ ^ ( mm m / r \ \ 
The left panel in Figure [4] shows the resulting structure functions of synchrotron intensity on the sky in this isotropic 




model for 7 = 2 and 7 = 4. While equal at large separations, two cases show small functional difi^erence in the 
transitional regime where effects of non-linear mapping between the magnetic field and synchrotron manifest them- 
selves. As the result, the small-scale amplitudes, where the synchrotron structure function essentially linearly reflects 
the mag netic field structure function, differ for 7 = 2 and 7 = 4 by the constant factor ~ 6/5, in accordance with 
Eq. (22). Such constant factor, however, is not readily observable if synchrotron correlations are mapped over small 
scales only. It is essential to track the behaviour beyond the scaling regime through non-linear transition to large 
scales. This encounters the following issue - the nonlinear, 7 sensitive, response of the synchrotron takes place where 
d cx / D{oo) ~ 1, i.e close to energy injection scale r/. Thus, to map it with accuracy sufficient to separate the 
different 7 cases, would require an accurate modeling of the largest scales of the turbulent cascade. 
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The right panel in Figure |4] demonstrates the power spectrum tha t correspond s to th e structure functions in the left 
panel. The power spectrum for 7 = 2 has been studied in detail in Chepurnov] ( 1998 ). The nonlinear mapping from 
the magnetic field to synchrotron emission affects the spectrum at long waves K rj < 10, proportional to the scale of 
turbulence saturation. The short scale asymptotics oc iir~^^+'"' follows the linearized relation between magnetic field 
and synchrotron statistics, with only the amplitude but not the slope being weakly dependent on 7. 

3.5. Correlation function of intensity for arbitrary 7 
The investigation in this section allows us to approximate the emissivity correlations by the ansatz 



ChT (r) « r{j)^H-^ (r), Dh-,, (r) « ^(7)7'(7)1?h^ (r) 



(35) 



where the strongly dependent on 7 amplitude 7^(7) 



HI 



is factorized from the scaling and angular 



dependences described by 7 = 2 term. For the structure function, adjusting the amplitude in a weakly 7-dependent 
fashion with ^^(7) ~ 1 makes the approximation even more accurate at small scales. 
This ansatz translates to the observable structure function of the synchrotron radiation 

DsyncjiR) « -4(7)7' (7)i?.y„c,7=2(R) (36) 

according to which the scaling and angular dependence of the synchrotron structure function can be un derstood from 
studying 7 = 2 case. For isotropic magnetic fields the above expression can be rewritten using Eq. ( 18 1 



DsyncM'^) ~ Aij)2''-^H^-'-* (r [1 + 7] - r [1 + ^7] ') 



.(R) 



(37) 



In case of axisymmctric turbulence one can choose the x coordinate to be aligned with the sky projection of the 
symmetry axis. In this frame covariance of the components of the magnetic field is diagonal, tr^j, = and the ratio of 
the variances ^(7) can be expressed via — a^x + (^yy and e = to give 



D,yr,,^^ (R) ~A(a)a'^^-^D,yr,c^^=2 (R) 



(38) 



1-e 



1 



(l-e2)2r(l+7)2Fi 



-, 1 



7 -, 2 



r 1 



2F1 



-, 1 



7 -, 2 



Eq. 



(|39| relates the structure function of synchrotron intensity for arbitrary index of relativistic electrons, i.e. for 

The additional 



arbitrary 7, with the structure function for synchrotron intensity fluctuations corresponding to 7 = 2 
factors do not depend on the distance between points for which the correlation is thought, but uniformly change the 
amplitude of the structure function. Both 7 and the e can be obt aine d independently as (see Discussion). Further we 
simplify out treatment focusing on the case of 7 = 2, while Eq. (39) allows us to generalize the results obtained for 
7 = 2 for an arbitrary index 7. 

4. ANISOTROPY OF THE SYNCHROTRON INTENSITY STATISTICS 

This section is focused on describing synchrotron fluctuations arising from a general model of axisymmctric 
anisotropic magnetic turbulence'^. On very general grounds one can expect mean magnetic field to influence statistics 
of magnetic turbulence. 

4.1. Power spectrum of the axisymmetric MHD turbulence 

The axisymmetric tensors presented in §2 describe structure functions of turbulence. In terms of the spectral 
representation of axisymmetric turbulence a general form of the index-symmetric correlati on tensor of the solen oidal 
vector field in the presence of a preferred direction A is given by two scalar power spectra ( Oughton et al.|[l997 ) 



(Hi Hi 



1 



{2n) 



E 



(k) (6 



F(k) 



(k • Xykikj + Xi\j - (k • X){kiXj + kjXi) 



l-(k-A)2 



(39) 



In the isotropic limit only E spectrum is present and is function only of the wave number, i.e. E{k). For axisymmetric 
turbulence, -E(k) and i^(k) functions depend on the wave number and the angle = k • A between the wave vector 
and the symmetry axis. The form of Eq. ( 39 ) is fully general if the statistics has a mirror symmetry in addition to the 
axial one ^. The normalization of tensors is chosen such as to have angle- independent traces. 

What combinations of E and F spectra correspond to independent modes is determined by the physical mechanisms 
of wave excitation and propagation. In the context of MHD, where the axis of symmetry is associated with the direction 

^ A particular model of MHD turbulence that corresponds to numerical simulations is consi dered in th e next section. 

* In the most general form two other pseudo-scalar spectra, C and H in the notation of Oug hton et al.| |l997), that describe correlations 
between toroidal and poloidal components of the magnetic field appear. Both terms are traceless and do not contribute to the rms energy 
of the magnetic field and are absent for axisymmetric turbulence with mirror symmetry. One of them, H, describes helical correlation 
resulting in index antisymmetric part of the correlation tensor which is irrelevant for synchrotron studies. We shall not consider possible 
contribution of the remaining C term in this paper. 
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(perhaps local) of the mean magnetic field, correlations due to Alfvenic modes correspond to the linear combination 
of "iJ-type" and "J "-type" spectra with F (k) = —E{k), while the slow and fast modes are described by the "F-type" 
correlation tensors ( Yan &: Lazarian ||2004', see also below). 



When anisotropy is global, the scaling of power with wavenumber k and the angular dependence on factorizes. 
Limiting ourselves to the power-law scaling we shall consider 

E{k,fi) ^ AEk-^'^^^Eifi) , F{k,^x)^ Apk-^'^'^Fi^i) (40) 

where m = 2/3 would correspond to Kolmogorov scaling. For example, in Alfvenic turbulence the fluctuations of the 
field are suppressed in the direction of the mean magnetic field, so one expects E (/i) to disfavor /i = ±1 direction. We 
normalize the angular functions to have unit monopoles in Legendre expansion, E{^) = ^'^q EiPi(^), etc., 

Eo = I J ^d^iE ifi) = 1, Fo = ^ J ^dfiF ifi) ^ 1 (41) 

In the presence of anisotropy the magnetic field components parallel and perpendicular to the symmetry axis have 
unequal variances, ajj = ((H • A)^) ^ ^a'j_ = ^((H x A)^). Since ((H • A)(H x A)) = due to azimuthal symmetry, crjj 
and |(T^ are the eigenvalues of the variance tensor (Jij — {Hi{x.)Hj{x)) , and are given by 

'yf\=^ f ^ (1 - M')^(/i) + ^ £^ df, (1 - ^,')F{^,) (42) 
= ^ £^ rfM (1 + f^'mn) + ^ £^ rfA* ^^'F{n) (43) 

where Pe,f = 4:TtAe,f J k'^dkk^^^"^'^-'' are measures of the power in, respectively, E and F parts of the spectrum. 
Note that individual amplitudes can, in principle, be negative, while the conditions Pe + Pp > and 4Pe + Pf > 
must hold. 

The ratio of the two variances is constant in the case of measurement performed in the global system of reference 



2^2 Pe + Pf-II PeE2 + PfK) 



Pe + ^Pf + Jo [PeE2 + PfF2 



(44) 



The parameter = 1 in the isotropic case when Pe = and E{fi) = 1. Values To- < 1 indicate the suppression while 
Ta > 1 the enhancement of the variance in the component parallel to the symmetry axis. Notable anisotropic limit 
of interest is that of Alfvenic turbulence, where Pf = —Pe and the component of the magnetic field along the mean 
is not perturbed, hence = and = 0. i?-type perturbations with completely suppressed waves along the mean 

magnetic field, Pf = 0, i?(/i) c>c (5(/i) lead to = 2. On the other hand, the essentially anisotropic F-type modes give 
To- = 4 even when the spectral function F is angle-independent. 

4.2. Transition from local to global axisymmetry 
In many physical situations, the preferred direction is local in nature and is only approximately maintained over 
the extended regions. In particular, in MHD turbulence A is associated with the direction of the mean magnetic field 
which may and usually does change over large scales. In this paper we are interested in the possible anisotropy of the 
synchrotron signal due to existence of a preferred direction over extended regions of space, i.e a global effect. Such 
anisotropy is expected when the direction of the magnetic field has a long-range coherence, i.e A can be thought of as 
fluctuating around the mean global direction Aq. 



One can estimate the effect of A "wandering" by averaging the tensor form of the correlations given by Eq. ( 39 ) over 
the distribution of A around the mean Aq (see wandering induced by Alfvenic turbulence in Lazarian & Vishniac 2000). 
Qualitatively the outcome is clear if the distribution of polar angle around Aq is isotropic. Then, at every level of the 
variance of A, the axisymmetric nature of the statistics is maintained, so the correlation tensor must still possess the 



form (39 1 with global direction Aq defining the symmetry axis. What changes when we average over local A directions 
is the distribution of power between E and F components, with power shifted to isotropic E term as the variance of 
A increases. 



(45) 



{H,ik)H*{k)) = E{k, k • Ao) [6,, - hk 
+F{k,k-\o) 



w n \ six ; ; ^ 1 w \ ^o) hkj + AqjAoj - (k • Ao)(fciAo-, + kjXoi) 

Wi{k ■ Ao, o-a) 5,j - kikj + I^l k • Aq, ax)- ' \ - 

j\ J J I l-(k-Ao)2 



where overline indicates the average over the cosine x = A • Aq. Wi and Wl are weight functions that depend on the 
measure of the average A wandering a\ so that WL(k • Aq, 0) = 1 and Wiik ■ Ao, 0) — WLik ■ Ao, oo) 
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While details of the averaging over local directions of A depend on exact form of the power spectra, main effect can 
be demonstrated by just averaging the basic tensor forms. In particular, the averaged i^-tensor (assigning 1 — (k • A)^ 
norm to the spectral function for convenience) becomes 



(k • + - (k • \){ki\j + kjXi) 



(46) 



r 



1) ((k- Xo)^kikj + XqiXoj - (k- Ao)(fciAoj +kj\Qi) 



We see that wandering of the local A, i.e. 7^ 1> leads to the appearance of the isotropic tensor part and the 
suppression of the original i^-type part. Although the exact weight functions between the parts depend on details, 
the structure of the transition described by Eq. 46 is general ~ the original tensor structure (of any form) that one 
starts with when A = Aq disappears as A gets distributed more and more isotropically, — J" 5 1 being replaced by the 
isotropic tensor. 

To estimate the characteristic angle of the transition to isotropic tensor structure we plot in Figure [5] our 
model weight functions Wj{(7\) = 5 (^1 ^ X'^^ a-nd WL{a\) = ^ ^3%^ — 1^ assuming a simple distribution P{x) = 

— s^P d,x- With this parameterization, ax matches very closely the mean deviation of cosine x from 

unity as long as a\ < 0.4, We find that one can approximately use the local tensor structure of the correlation for 




COS 

n/S n/4 37r/S n/2 

Fig. 5. — Example of the transition weights between anisotropic and isotropic correlation tensor structures as functions of the mean 
angle between a local and the global symmetry directions, cos~^ x = cos~^ ( A ■ Aq ) • 



If the mean fiuctuation of A exceeds « Stt/S, one 



< cos ^ X ^ 37r/8 one has a mix of the local and isotropic 



global studies when the mean deviation of A from Aq is less than tt 
can use the isotropic tensor. In the intermediate range tt 

contributions. Two contr ibut ions are in similar proportions when cos~^ x ~ 7r/4 

The averaging in Eq. (45) is an ensemble averaging for three dimensional correlations. It worth noting that the 
averaging arising from field wandering in Eq. (45) acts differently on the different parts of the correlation tensor. 



4.3. Structure Junctions at small scales 

We have argued that the correlation properties of the synchrotron emission at small scales are well approximated 
by the linearized limit of the 7 = 2 model. This allows one to study scaling properties of emission for 7 = 2 and 
then generalize those for any astrophysically relevant 7. The full expression for the structure functional for 7 = 2 (see 
Appendix ICl) 



rfzii/i(Xi,zi) 



dz2Hl (X2, 2:2) 



(47) 



dz+ J dz {4a,, [Dj, (X, z) - Dj, (0, z)] - [D,, (X, z) Aj [X, z) ~ Dj, (0, z) D,, (0, z)]} . 

in the small-scale limit (i.e. when the structure functions are much less than dispersions) provides a linearized expression 

DsyncCR) - 4ct„- J dz [Dj, (R, z) - Dj, (0, z)] . (48) 

that involves x, y — 1, 2) components of the magnetic field correlation tensor. 

For the following discussion we introduce the coordinate frame that has z-axis along the line-of-sight and y-axis 
perpendicular to the plane spanned by the line-of-sight and the symmetry axis A. The x-axis is then along the 
projection of A onto the image {x,y) plane, denoted by the unit 2D vector A. The vector A now has components 
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Ai = (sin(?,0,cos^^). The polar angle of the 2D separation R is denoted by </> while 'tp is the polar angle of the 2D wave 

vector K. Both have the meaning of angles between the correspondent vectors and A. 
The 2D part of the dispersion tensor dij is diagonal in the chosen frame 

o-xx = o-y sin^ d + cos^ , (Jyy = , a^y = . (49) 
Its anisotropy can be represented by the parameter 

^^ <T,.-<Tyy ^ sin^g(af-iai) ^ sin^g(r.-l) 
cTa:x+(Tyy sin^ ^ (crf - ic72 ) + smH{r^-l) + 2 

which depends both on the properties of 3D turbulence and the projection angle 9. The parameter e play a key role 

in interpretation of the synchrotron signal. 

The integration over the line of sight is equivalent to setting fc^ = in the spectral domain. The z-axis projected 
correlation function of the field becomes 



Jdz{HiHj} oc Jd^Ke''^-^ E{K) (jy - A', A',) +sin^9 F(K) 



cos^ipKiKj + AjAj - cos'tp{KiAj + KjAi) 
1 — cos^Vsin^^ 



(51) 



where, as we discussed earlier -E(K) and -F(K) functions describe the underlying statistics of turbulence. 

The 2D power spectra E{K) ^ E{kz = 0) and F(K) = F{kz = 0)^ depend, parametrically, on the angle 9 between 
the symmetry axis and the line of sight. Notably, the contribution from the F-modes vanishes if the axis of symmetry 
is along the line-of-sight, sin^ = 0. 

More interesting is the dependence on the positional angle tp. One can deal with the harmonic decomposition of 2D 
spectra which is defined as En{K) ~ ^ Jo'^dip e^"^^ E{K,tp), and analogously for Fn{K). 

Below we obtain expressions for multipoles of the decomposition of 2D structure functions of magnetic field. These 
multipoles can be studied observationally. In what follows we transform back to the correlation function domain we 
express the spectral tensors by the differential operators. We consider separately the additive E and F contributions 
to the structure functions. 

4.3.1. Multipole expansion for E-term 
Correlation function that arises from S-type spectral tensor can be expressed as 



^ dz{HiHj) oc {didj - SijA) $B(i?, (/.) , (52) 
with 

^e{RA)= Y1 i'^e''"^ K-^dKJn{KR)En{K) (53) 

n— — C30 

The synchrotron structure function has a similar form 

DsynciR) oc 8(7,j {d^dj - (5y A) ^reg{R, 4>) (54) 

but with the regularized ^reg{R, 4>) = <I'(0, 0) - ^(i?, 0) It becomes 

Dsync{R, (t>) -(o-xx + (yyv)A^reg{R, <P) + {cTxx - (^yy) {d^ - dy) '^>reg{R, <P) ■ (55) 

Both terms have concise multipole expansions 



A^reg{R,cl>)= ^ iV"M AT^iAT [J„(i^i?) - J„(0)] 
n=-^oo -^0 

oo ^ oo rp -I- P 

{dl - d^y) <^reg{R, <^) = E ^"^""^ / ^'^^ lUKR) ' J„(0)] ^ 

that lead to the expression for the multipoles of the normalized synchrotron structure function 

D^{R) = ^ J^\^e-'"^D,y^,{R,^) = 

= r / KdK [J„(0) - Jn{KR)] [En - |e (£„+2 + En-2)] ■ 
Jo 



(56) 
(57) 



(58) 



® We shall use the same letter notation for 3D and 2D spectra, which, hopefully, should not lead to the confusion. 

® We have left 4> dependence in ^>(0, (j>) to signify that regularization may aifect not only the monopole but also the higher harmonics. 
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The first term directly maps the multipole terms of the power spectrum to correspondent multipole terms in the 
projected correlation function, while the second one describes the multipole coupling in the presence of anisotropy. In 
particular, the spectral quadrupole feeds into the monopole of the structure function 

/•oo 

Do{R)^ KdK[l-MKR)][Eo-le{E2 + E^2)] , (59) 
Jo 

while the quadrupole dependence of the correlation is affected by the monopole and the octupole of the spectrum 

D±2{R) = / KdKJ2{KR) [E±2 - (Eo + E±i)] . (60) 
Jo 



With factorization Eq. (40) in place, the multipole terms of the synchrotron intensity in small scale regime and 7 = 2 
are 



D.r, {R)~AECn{m) (e„ - \e (^„+2 + ^„-2 ) ) , (61) 



The expression for the coefficients is valid for m < 1. For the Kolmogorov index, m — 2/3, Co(2/3) ~ 1.12, 6*2(2/3) sa 
0.51, 6*4(2/3) « —0.03, decreasing to under 1% of the monopole value for higher n. Lower m somewhat enhances the 
higher muhipoles, e.g. for m = 1/2 Co(l/2) w 0.93, 6*2(1/2) w 0.40, 6*4(1/2) w -0.04, although they remain smafi 
beyond n — 4. 



Eq. ( 62 1 relates the multipoles of correlation function that can be measured to the multipoles of the underlying 



power spectrum that characterizes axisymmetric turbulence. The dependence on the spectral index of turbulence 



enters both through _R^+'" and through 6*„(to) terms given by Eq. (62) 



4.3.2. Multipole expansion for F-term 

Expression for the F-type correlation function is somewhat more complex due to explicit angular dependence in 
the spectral tensor. While it is possible to represent it as the fourth-order derivatives of a single scalar function (see 



[Oughton e t alj 19971 here we shall use the direct coordinate approach. 
From Eq. (^5lL two relevant component correlations are 



1 — cos^f/j siia^O 
sin^O f ,2r^ iK R T^/r^\ cos^V siu^ V 



J dz{HyHy)^'^ J d'Ke^^-^F{K) - 



cos^V' sin^ 



Combining these expressions into Eq. (48) for DgynciR, 4') £^nd performing multipole expansion we fi nd that the 



multipole coefficients of the synchrotron structure function from waves of i^-type are (compare with Eq. (58l) 



/>oo 

D^{R) ^ i" sin^ 6 KdK[Jn{0)~Jn{KR)]J2[Pp-hiPp-2 + Pp+^)]Gn-p (63) 
where harmonic coefficients of the angular function that describes the trace, xx + yy, of i^-tensor (divided by sin^ 6) is 



27r 7o 1 — cos^ip sin^ 

This function describes how the spectrum contributes to different angular multipoles in the observed synchrotron 
correlations. It is especially important in many interesting cases when the spectrum F{K.) is isotropic, since then it 
determines the angular behaviour of Dsync in its entirety. 

Geometric functions that arise in our studies are presented, in particular graphically, in Appendix |D] In two limiting 
cases, contains only low harmonics. For 6* = ^, 6*^ contributes the monopole 6*^(0) — ^ and the quadrupole 
G^2'(0) = ^5^0^(0) terms only, while in the opposite limit, 6 — ^tt, there is just the monopole Gq{tt/2) = 26*^(0). 
Varying the symmetry orientation angle 6 (see Figure in Appendix [d]) changes the power distribution between har- 
monics Gp (9) in a rather complex way due to singular behaviour a.t 9 ^tt. At steep angles < 7r/3 as symmetry axis 
approaches sky direction, the whole train of high multipoles is generated, albeit at decreasing amplitudes. This may 
or may not lead to observable consequences, especially in the context of steep spectra that suppress high multipole 
contributions even further. 

The relation between multipoles of the observed intensity correlation function and multipoles of the underlying 
spectrum in case of the _F-term is 

00 

D„{R) ^ AFCnim) sin^e ^ [Pp - (Fp_2 + Fp+2)] G^-pi?'+™ • (65) 



We reiterate that if the symmetry axis is along the line-of-sight this signal is not observed due to amplitude suppression from linc-of-sight 
integration. However, this case is relevant for consideration of the fast modes in high-/? MHD turbulence 
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5. SYNCHROTRON EMISSION FROM MHD TURBULENCE 

In the previous section we dealt with a general form of the tensor of the axisymmetric solenoidal field. A realistic 
MHD turbulence corresponds to the particular form of this tensor. In what follows we use the model of compressible 
MHD turbulence that follows from numerical simulations and also can be justified theoretically. Predictions obtained 
within this model can be used to test the model and can provide insight into the role of compressible (fast and slow 
modes) versus non-compressible (Alfven modes) motions. 



5.1. Model of Alfvenic turbulence 

It has been well known that small amplitude magnetic perturbations can be decomposed into Alfven, slow and fast 
basic MHD modes. However, as the conventional discussion involves the decomposition in respect to the mean field, 
such a decompositions becomes questionable for any magnetic field pertu rbations of appreciable amplitude. The use 
of local system of reference (see Lazarian & Vishniac 1999 Cho &: Vishniac,,2000J makes the p r ocedure justifiable. 
The s tatistical decomposition of the modes has been performed in Cho & Lazarian ( 2002 2003 ) ; Kowal & Lazarian 
(2010). It was proven that the modes provide cascades of their own and the drainage of energy between the cascades 
is limited. 

Alfvenic modes are the most importan t ones with the mo st of the energy residing in them and they are less subject 
to damping in astrophysical plasmas (see Cho fc Lazarianj2002). They have the statistical properties which are similar 
to the Alfvenic modes obtained in incompressible simulations ( Cho fc Lazariaii||2003 ). This motivates us to study the 
Alfvenic turbulence as our primary example 



The spectrum of Alfvenic modes was obtained in Cho et al. ( 2002 ) for the driving turbulence at large scale with 
equal to the Alfven velocity: 

ku 



S(k)0C/fc-l"/3i-l/3exp 



2/3 
± J 



(66) 



where L is the injection scale of the turbulence. The exponent in Eq. (66 1 reflects the so-called critical balance condition 
between the parallel and perpendicular modes. It shows that most of the energy resides around the wavenunibers where 
the turnover of the turbulent eddies is equal to the period of the Alfven wave. 

Note that the driving at Vl ~ Va corresponds to model of GS95 turbulence. In the case of weaker driving at the 
outer scale, the turbulence becomes sub- Alfvenic at all scales and the relations between fcy and k± from LV99 should 
be used. Namely, the critical balance condition provides a modified relation between the parallel and perpendicular 
modes at sufficiently small scales 

fc|| ^L-i(fcii)2/3M^/3 (67) 

where Ma is the Alfven Mach number, which is Vl/Va- Therefore one expects that the description of the turbulence 
for Ma < 1 involves the ratio of the left and right sides of the Eq. ( 67 ) in the exponent of Eq. ( 66 1 . This would provide 
the description of turbulence in the local system of reference. 

The tensor structure of Alfven power spectrum is derived from the condition that displacement and velocity of 
plasma in an Alfven wave k are orthogonal to the plane spanned by the magnetic field and k, va cx k x A. Magnetic 
field response is given by the condition that fluctuating field is frozen. 



H{k) (xkx [v X Xj /uj{k) , 

where for the Alfvenic waves the frequency cj oc k • A . Simple algebraic manipulations then show 

(k • X)'^kikj + XiXj - (k • X)(kiXj + kj\) 



(i/.(k)H;(k))^« (4 



i-{k-xy 



(69) 



which in terms of the general E - F decomposition (391 corresponds to Fik) = —Eik). This tensor refiects the absence 



of the perturbations in the magnetic field component parallel to the mean field A. 

In the global system of reference, however, the anisotropy is scale independent and determined by the anisotropy at 
the outer scale. In this system of reference the tensor of magnetic perturbations can be given by 



E{k) cx fc-"/3 exp 



k-Ao 



(l-(k-Ao)^ 



2/3 



(70) 



arises from extending to global scales and Ma 7^ 1 the phenomenological, locally anisotropic, model of Cho et al. 
(2002). The axis of symmetry is associated with the direction of the mean magnetic field Ao. The spectrum reflects the 
consideration that in critically balanced turbulence with Ma < 1 the modes with parallel to magnetic field wavenumbcr 

fc|| that exceed the critical value ~ L~^/'^k'^'^ Ma^^^ are suppressed. At the injection scale L, if fc|| '-^ L. 
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Following Section 4.2 the global tensor structure is a weighted mix of the local tensor of Eq. 69 and the isotropic 
tensor form. The — section of the spectrum becomes 



(F,(K)ff;(K)) cxi^-"/=^exp 



((5,, - k,k,^ Wi 



|cos^/'| smO 




KiKj I — sin 



1 — cos^ ip sin^ 9) "^^^ 

2^ cos' tpkikj + AjAj - cos'4!{KiAj + Kjk^i 



1 — cos^^j sin 



Wl 



(71) 



Now A is the projection of the global symmetry axis Aq but for brevity we are omitting index '0' henceforth. We 
note that if the symmetry axis is strictly along the line of sight, sin 6* = 0, the anisotropy effects are not observable 
in sky projection. Maximal anisotropic effects are observed when the symmetry axis is perpendicular to the line of 
sight, sin0 = 1. At small Mach numbers, parallel to the mean magnetic field modes are strongly suppressed and the 
turbulence is hig hly anisotropic. The effect is retained in sky projection, but is mitigated by sin0. 

The model Eq ( 71 ) that we use to discuss Alfvenic turbulence has two parameters - the Alfvenic Mach number Ma 
and the mean deviation angle x that gives the weights W^l(x) and W/(x)- On the other hand, the level of magnetic 
field wandering due to turbulence is itself related to the Alfvenic Mach number. High Ma > 1 turbulence leads to 
significant bending of magnetic field lines and isotropization of the observable statistics. Conversely, details of the 
local mode structure of the correlation tensor can be observable only in sub- Alfvenic, Ma < 1 regime. As the guiding 
model of the wandering consistent with given Ma we shall employ a rough prescription that Ma gives a tangent of 
the mean deviation of the local angle from the global axis, x ~ 

1/^1 + ^2. The sub- Alfvenic regime but with some 
isotropization may be relevant to synchrotron studies if the line-of-sight integration averages over regions with different 
magnetic field directions. 

Asymmetry in the variance in sky components of the magnetic field in this model is shown in Figure|6] The parameter 
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Fig. 6. — Asymmetry in the variance of the Hx (parall el t o sky projection of the mean field) and perpendicular Hy, components of the 
magnetic field for the model of Alfvenic turbulence, Eq. | |70[ l as parameterized by e(9) = (uxx — (^yy) I {pxx -\- i^yy), as the function of the 
projection angle 9. Left: change from the sub-Alfvenic (Ma = 0.1 bottom solid) to trans-Alfvenic (Ma ~ 1, top solid) to super-Alfvenic 
{Ma > 1, middle solid) behaviour. Averaging over local direction of the magnetic field is carried out in accordance with Ma- Right: 
dependence on Mach number Ma when field wandering is neglected. The bottom solid curve in the negative sector is Mach independent 
result (ro- = 0) when local tensor correlation form is assumed to hold exactly, e > sector shows the dependence on Mach number 
if isotropic i?-tensor (but with anisotropic E(k.)) is used to model Alfvenic correlations. The dashed curves, in decreasing order of e, 
correspond to Ma = 0.2, 0.5, 2. The limiting behaviour with the highest anisotropy e ( top solid curve marked Ma <SC 1) is closely followed 
already at Ma < 0.2. 



e, given by Eq. (50), as expected, is larger by magnitude when the symmetry axis is aligned with the sky plane, and 
is vanishing when it is parallel to the line of sight. Left panel shows that e magnitude is small to moderate when 
field wondering is appropriately accounted for. e is slightly negative for strongly sub-Alfvenic turbulence, almost zero 
at Ma ~ 0.4, reaches maximum positive values, e ~ 0.1 in trans-Alfvenic regime M ^ 1 — 2 and then decreases to 
isotropic limit e — >■ for super Alfvenic Ma- However the result is somewhat sensitive on exact procedure we model 
the transition. 



The right panel provides a cautionary illustration. If one adopts the exact local description Eq. ( 69 1 for the Alfvenic 
turbulence to hold through to global scales, Wj ~ 0, Wl « 1, then we get negative and practically~Ma-independent e, 
since ^ 0. This means that the perturbations in magnetic field component parallel to the mean field are suppressed. 
On the other hand, if anisotropic power spectrum is used with isotropized £'-tensor form to describe the turbulence 
the result is opposite - parallel component is presumed to fluctuate more than the perpendicular one and < e < 1/3, 
with increasing Ma taking us to isotropic limit. Far from being just of academic interest, this latter case appear in 
the description of a mix of Alfven and Slow modes in incompressible limit. We also see that averaging of the statistics 
over local field direction may effectively bring us to this regime, but probably not at the highest magnitude of the 
effect given by = 2 and e w sin^ 9/ (sin^ 9 + 2). 

In Figure [T] the distribution of power in multipoles i?„ is shown for representative Alfvenic Mach numbers and several 
orientations 9. The multipole expansion contains only even multipoles and the real nature of coefficients reflects the 
presence of only the cosine terms in the decomposition. The anisotropy of the spectrum is predominately a quadrupole 
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Fig. 7. — Multipole distribution of power in the units of the monopole En/ Ea for Ma 
(circles, squares and diamonds, respectively). 



1 (left) and Ma = 0.4 (right) for I 



one for Ma > 1. As Ma decreases, the higher multipoles appear and in the hmit the spectrum contains all En of equal 
amplitudes with compensating alternating signs. The presence of high m ultipoles is largely due to sharp features in 
angular dependence of the spectrum which are artifacts of the model (see Cho et al. |2002, for technical remedies of 
this issue) 



Th e st ructure function o f th e synchrotron fluctuations from Alfvenic modes is a special combination of Eqs. (611 
and ( 65 ) that matches Eq. ( 71 1 



Wi 



En — 2^1 En+2 



En- 



Wl 



E 



E.^ 



p-2 



'p+2 



G 



71 — p 



R^/^ . (72) 



where we refer to Appendix [d| for the properties of G^{9) function. In many cases Alfvenic modes are expected to 
dominate magnetic field perturbations and therefore Eq. ( 72 1 provides a good representation of the multipoles of the 
structure function of the synchrotron intensity for small 2D separations between the lines of sight \R\ <^ L, where L 
is the injection scale of the turbulence^. 

5.2. Fast modes 

Fast modes are compressible modes which in high beta plasma, i.e. plasma with gaseous pressure much larger than 
the magnetic one, transfer to sound waves, while in low beta plasma propagate with Alfven velocity due to compressions 
of the magnetic field. These modes marginally interact wit h Alfven modes and were show n to form their own cascade. 
They were identified with an isotr opic spectral function in Cho & Lazarian (2002, 2003). The tensor form of the fast 
mode turbulence was presented in Yan & Lazarian ( 2004 ) 



k-y^ (k • X)^kikj + X^Xj - (k • X){kXj + kjXi) 



F 87rLi/2 



1 (low P) 

1 - (k • Xf (high P) 



(73) 



1 - (k- A)2 

The difference between fast modes in high-/3 and low-/3 plasmas comes from different behaviour of the velocities in the 
medium. While by very definition of the fast modes the velocity at a wave k lies in the plane spanned by k and A, 
for /3 ^ 1 velocities are potential, while for /3 <C 1 they are orthogonal to the direction of the magnetic field A. The 
frequency of fast modes Wi?(k) has a weak additional dependence on the direction of the wave vector if /3 ~ 1, the 
effect we do not consider here. 

In contrast to Alfven modes, fast modes are of F-type. Despite the isotropic power scaling, fast modes are not 
statistically isotropic, with anisotropy built into the tensor structure of the spectrum. 

Asymmetry in the variance of the sky components of the magnetic field for fast modes is plotted in Figure [Sj 
Compared with that of the Alfven modes, fast modes always produce positive anisotropy e. This effect is quite distinct 
from Alfvenic modes where e is much smaller, or even negative. Determining e from fitting multipole composition of 
synchrotron correlations may give result that distinguishes whether we are dealing with fast (positive e > 1/3) modes. 



Isotropic nature of the power scaling leads to a very simple application of Eqs. (651 and (61) to MHD fast modes 



In low /3 case we have -F'(K) oc K which corresponds to m = 1/2 and only the monopole Fo present. This gives 
bn{R) - ApCniim [Wi {5nQ " eJ„2) + Wl sin^O {G^{9) - \e {Gl_^{d) + G^+2(^)])] % R^''' {low /3) (74) 

Additional angular dependence in the spectrum in high-/? medium can be accounted for by just evaluating the 
coefficients at 6* = 

Dn{R) ^ ApCnil/i) [Wi {Sno - eSn2) + Wl sin^e (G,^(0) - [G^_2(0) + G,^+2(0)])] Fo R""^ {high p) (75) 

® For the observed correlation functions this means that the angle <l> <C L/ D, where D is the distance to the object. 
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Fig. 8. — Left: e{6) = {ctxx — o-yy) / {f^xx + f^yy) I as the function of the projection angle 6 for the fast modes. Solid lines: change in 
high-/3 plasma from the local behaviour Wj^ = 1, Wj = (top) to the isotropized limit = (bottom) through intermediate cases that 
corresponds to average angle of local A deviations of 7r/8, 7r/4, 3n/8 (from top to bottom). Dashed curves - same for low-/? plasma. Shading 
covers the interval that different /3 will fill for the same cos^^ x- Right: same for Slow modes. Isotropization of the slow modes is treated 
as.associated with the increase. in the Alfvenic Mach number., ... , . „ , . i. • i i r 

We note xnat m case of nign p Ine averaging over A wandering is exact, so it is useiul to write explicitly the set ot 
all non-vanishing multipoles the synchrotron structure has in this limit. Taking into account the signs of C„(l/2) 
coefficients we have 
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(76) 



The structure of the multipoles of correlation functions is different from that of multipoles arising from Alfven modes, 
which allows one to separate them, finding the measure of media compressibility. 

5.3. Slow modes 

Slow modes are usually subdominant but provide an interesting case of the -F-type modes but with the power 
spectrum that follows the Alfvenic power spectrum. Displacements in slow modes lie in k - A plane and are orthogonal 
to the Fast mode ones. In the low /3 limit, the slow displacements are parallel to the magnetic field and do not lead 
to magnetic field perturbations. So we focus at /3 ^ 1 case where the slow mode motion of plasma is perpendicular 

to the wave vector, w oc k x ^k x A^ The frequency of the slow mode in this regime is similar to that of Alfven ones, 
ujs cx k • A. The frozen field condition gives the magnetic field correlation tensor of the F-type 

(k • Xfkikj + Xi'Xj - (k • X){hXj + k^Xi) 



(Fi(k)F;(k))^ ^ i?(k) 



(high /3 only). 



(77) 



l-(k-A)2 

Figure |8] shows that the anisotropy in the sky component of the magnetic field variance produced by Slow modes 
is similar to the Fast modes demonstrating isotropization as Ma increases. We stress, however, that e is determined 
by all the modes excited in the medium and cannot be simply separated in individual contributions. The e xpre ssion 
for the synchrotron structure function multipoles can be found by combination of the previous results Eq ( 65 ) and 
Eq. ffol) 
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With slow modes being subdominant it is more relevant to study their contribution to synchrotron in combination 
with other modes. The largest effect the slow modes have are in /3 ^ 1 limit that describes, in effect, an incompressible 
plasma. Slow and Alfven modes in this limit are two "linear polarizations" of a transverse displacemen t wave. In 
particular, if in strong turbulence the powers in slow and Alfven modes are equal ( Cho fc Lazarian||2003 1 , the waves 
are unpolarized. The correlation tensor that describes such waves is of a pure i?-type 



(i/,(k)i/;(k)) 
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(79) 



In Appendix [B] we discuss the multipole structure of the 3D corr elation function that corresponds to such spectra. 
The synchrotron correlation structure that corresponds to Eq. ( 79 ) has the following multipoles 



which presents the structure different both from multipoles of Alfven and fast modes. 



(80) 
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5.4. Angular structure of synchrotron correlations from MHD modes 

Full observed synchrotron signal reflects the combination of Alfven, Fast and Slow modes in the proportion de- 
termined by the distribution of power b etw een th e m odes. T he s tructure function of the resulting emission can be 
obtained by linear combination of Eq. (72), Eq. (75) or Eq. (74), and Eq. (78). However the dependence on the 
underlying models is non-linear, since the anisotropy of magnetic held variance e is determined by all the contributed 
modes at the same time. 

First we observe that the Fast modes lead to different scaling of the synchrotron correlations than the Alfven or Slow 
modes, oc R^^^ versus cx i?^/'^, respectively. Here we suggest to focus on the measurement of the anisotropy of D{R, (j)) 
as independent from scaling signature of the contributions from different modes. Measurement of the moments of 
Dn{R) / Dq{R) may prove more robust discriminant between the contribution modes in cases when power-law scaling 
is either not exact or its accuracy is affected by observational algorithms or theoretical approximations. 



We first investigate the limit when only one of spectral contribution is dominant, either Alfven (Eq. 72), or Fast 
(Eqs 75][74 ), modes, thus assur aing e determined by the dominant mode only. Then, as an example of a mix between 
modes we discuss the case Eq. ( 79 1 of strong incompressible turbulence which maintains both Slow and Alfven modes. 

Figu re [9| shows the details of me quadrupole and octupole contributions for Alfvenic turbulence. Within our model 
of Eq. ( 72 ) the predicted quadrupole depends on the Alfvenic Mach number Ma and follows one of the curves in the 
left panelof Figure [9j 
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Fig. 9. — Left: Quadrupole of the synchrotron structure function in the units of the monopole versus the mean magnetic field to the 
line-of-sight angle d for the Alfven modes. Effect of Alfvenic Mach number, including expected isotropization due to fluctuations of the 
local field direction. The curves correspond to Ma = 0.1,0.4,0.7,1,2 ( from lower to upper ones). Right: same plot for the octupole. 



First observation from Figure [9] is that the expected quadrupole of the synchrotron correlations is negative. With 
structure functions scaling as ~ (Dq ~\- D2 cos (j> . . negative quadrupole means extended iso-correlation contours 

along the x-axis, the sky projected direction of the magnetic field. In case of Alfvenic turbulence this corr esp onds to 
spectral suppression of the modes parallel to the field having the dominant effect. It is easy to see from Eq. (72 1 that if 
the power distribution was 3D isotropic, the local tensor structure of the Alfven modes would have produceoa positwe 
quadrupole in synchrotron, i.e the contours extended orthogonally to the symmetry axis. Thus, observing negative 
synchrotron quadrupole in Alfven turbulence context is a strong argument for anisotropic underlying 3D cascade that 
suppresses parallel modes. 

For our calculations we have adopted the model description that takes into account the fluctuations of the local 
direction of the magnetic field at the level determined by the Ma- Thus, changing involved both the field wondering 
effect and the change in the power spectrum. As the result of field wondering the anisotropy of the observed synchrotron 
becomes suppressed when mean field is roughly aligned with the line of sight, 9 < 7r/6, unless the turbulence is strongly 
sub- Alfvenic. If field axis is closer to be orthogonal to the line of sight, the anisotropy is strong. Alfvenic turbulence 
leads to the negative quadrupole which magnitude is hmited to \D2/Do\ < C2(2/3)/Co(2/3) = 0.45 as Ma decreases, 
with values close to the limit \D2/Do\ ^ 0.4 already exhibited at Ma = 0.4. while ID2I/D0 ^ 0.3 for trans-Alfvenic 
Ma ^ 1. Measuring D2/D0 « 0.4 will be a strong indication of the sub Alfvenic regime. Although diminishing with 
higher Ma, quadrupole is still significant, \D2/Di)\ ^ 0.1 even for mildly super-Alfvenic Ma ^ few. 

Right panel in Figure |9] demonstrates that the octupole contribution is quite small, at 2% [04^/ Dq ^ 0.02) level 
under the most favorable parameters. Indeed the steep spectra that are produced in MHD turbulence, while reflecting 
the multipole decomposition of the spectrum, demonstrate suppression of the multipoles higher than the quadrupole. 
This is the case for both Alfvenic and Fast modes and represents a challenge to measure them. 

Fast mod es h ave isotropic power spectrum and produce the anisotropy solely due to the anisotropy built into F- 
tcnsor, Eq. (73), which is somewhat different in the cases of high and low /3 plasma. Predictions of the quadrupole 
from the Fast modes are given in the left panel of Figure [TO} There are two parameters, /3, and the mean deviation of 
the local field direction At each x, for which three are shown, cos~^ x = 7r/16, 7r/6, 7r/3, D2/D0 from Fast modes lie 
distinctly in a narrow band limited by /3 ^ 1 and (3 <^ 1 predictions according to Eq. ( 75 ) and Eq. (l74| respectively. 



The Fast mode quadrupole anisotropy is rather uniform as the function of the orientation angle at 9 > n/G and at its 
maximum is smaller in magnitude than that from Alfvenic modes. Fast mode signal becomes nearly isotropic as field 
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Fig. 10. — Quadrupolo of the synchrotron structure function in the units of the monopole versus the mean magnetic field to the line-of- 
sight angle 9. Left; Fast mode contribution shown in bands for cos~^ x = '''/16, n/6, 7r/3, D2/DQ, decreasing in magnitude with increased 
field wandering. Each band corresponds to the range of different plasma /9, bounded by /3 2> 1 and /3 <C 1 limits. Right: the same quantity 
for the mix of Alfven and Slow modes in the regime of strong incompressible turbulence. The curves correspond to Ma = 0.1, 0.4, 0.7, 1, 2 
( from lower to upper ones). 



wandering increases, being only at 5% when local wandering is similar to Ma = 2 of Alfvenic case. The synchrotron 
quadrupole from Fast modes is never expected to exceed \D2/ D{)\ w 0.25. Thus, measuring higher quadrupole will put 
the importance of the Fast modes under question. 

As the last exercise we are considering the mix of Alfven and Slow modes in the regime of strong incompressible 
turbulence. This is the regime of pure E — type modes with Alfven like power spectru m d istribution and 3D isotropic 
tensor spectral tensor as modeled in Eq. (79 1. The results in the right panel of Figure 10 shows that the synchrotron 
quadrupole structure in this model is very close to that of Alfvenic turbulence. This confirms again that exact local 
structure of the spectral tensor in Alfvenic turbulence is not decisive in determining the observable quadruple of the 
synchrotron. It is primarily reflects the anisotropic distribution of power in the cascade. 



6. DISCUSSION 

Our paper is aimed to provide the quantitative analytical description of the statistics of synchrotron fluctuations. 
In this section we discuss what our results mean in terms of observations and what is the relation of this work with 
previous related publications. 

6.1. Our advances in describing synchrotron intensity statistics 

Synchrotron fluctuations can be characterized statistically by the slope of their spectrum and their amplitude. In 
this paper we have obtained expressions for the synchrotron fluctuation amplitude and have shown that the spectral 
slope of fluctuations depends weakly on the spectral index of relativistic electrons. The latter presents a substantial 
advance compared to the earlier studies, which were providing expressions for just one electron spectral index, namely, 
7 = 2 and were not able to estimate the errors arising from this assumption. Our approximation for correlations at 
arbitrary 7, namely. 
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where = a^x + cfyy are magnetic fleld dispersions and e = is the degree of anisotropy (see more in §3, 

Eq. (81) ) allows to express the statistics of fluctuations through the statistics obtained for 7 = 2. This expression 
opens avenues for analysis of the synchrotron intensities for a variety of astrophysical circumstances, for which 7 
substantially deviates from 2. In other words, synchrotron intensity studies for which possible statistical description 
was earlier limited to a very special case of a single and unique spectral index of cosmic rays has been extended to 
a wide range of spectral indexes which encompass astrophysically important situations. Eq. (81) is our first major 
result. 



Earlier studies were limited mostly to hypothetical isotropic magnetic turbulence (see Chepurnov|1998" ). This would 
provide a very special and rather unrealistic form for Dsync,-y=2 in Eq. (81). Thus the second part of our paper was 
devoted to obtaining adequate forms for Dsvnc.-/=2- 



Unlike other earlier studies (see|Lazarian fc Shutenkov|1990 ), we avoided dividing magnetic field into the regular and 
turbulent components. Such a division is unrealistic in view of modern day understanding of turbulence. Instead, in the 
presence of mean magnetic field, we adopted the model of the axially symmetric magnetic turbulence, where the mean 
field is directed along the direction of the axis. Both theory and simulations are suggestive that this description should 



correspond wel l in the global frame of reference, related to the mean magnetic field (see Lazarian & Vishniac"1999 



Cho & Vishniac 2000; Cho et al. 2002 1. To describe the axially symmetric ma gnetic turbulence we used the genera 



expressions m Dughton et al. "(1997), which generalize the classical results in Batchelor (1946) and Chandrasekhar 
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(1950). With this, our expression provided a combination of two different tensors with two functions E and F in front 
of them. At this point one can estabUsh the existence of the anisotropy of synchrotron signal and its relation to the 
aforementioned functions. 

To move one step further, we adopted a model of compressible MHD turbulence from Cho & Lazarian (2002 20031 
and expressed the statistics of the Alfven, slow and fast modes in terms of the combinations of E and F functions, 'i'his 
allowed us to express the expected anisotropics of the synchrotron emission through the physically motivated basic 
modes existing in MHD turbulence. Our approach allows to predict the synchrotron emission statistics if the turbulence 
is defined and, alternatively, study fundamental properties of MHD turbulence analyzing spectra and anisotropics of 
the synchrotron statistics. 

While we believe that the decomposition into modes in Cho & Lazarian ( 2002 2003 1 is well motivated theoretically 
and proven numerically, our approach in the present paper can be applied to any arbitrary model of axisymmetric 
turbulence. For instan ce, it can be applied to test the predictions of the model of MHD turbulence consisting from slab 
and 2D fluctuations^ ( Matthaeus et al]|l990 Zank fc Matthaeusl|l993 ). We are far from assuming that ast rophysical 
magnetized compressibl e turbulence is completely understood, -bbr instance, recent numerical simulations in |Federrath| 



et al. (2008 2009 



2010 ) demonstrate the effects of turbulence injection on the resulting turbulence. In the presence of 



multiple sources of energy injection, instabilities and damping processes the spectrum of magnetic turbulence inferred 
from synchrotron observations may be rather complex. Our description accounts for turbulence compressibility and 
our analytical expressions can be used for turbulence with several injection scales. This opens avenues for realistic 
magnetic turbulence in its complexity. 



6.2. Studying magnetic field direction 

Our results show that magnetic field direction can be studied through the analysis of the anisotropy of synchrotron 
intensity correlation functions. This provides a synergetic way of studying magnetic fields. 

Other approaches to studying magnetic field direction are known in the literature. Some of them are based on 
direct measurements of synchrotron polarization (see Ginzburg 1981), polarization from aligned grains (see review 
by Lazarian 2007 and references therein), or aligned atoms (see Yan & Lazarian 2008). In addition, a statistical 
determination of the mean magnetic field direction is p ossible with the statistical analysis of the Doppler shifted lines 



(Lazarian et al 



2002 Esquivel & Lazarian 2005 2011), that has similarities with the one we proposed here. In this 



technique the difference in velocity correlations parallel and perpendicular to mean magnetic field also reveals the 
direction of ma gnetic field. The a forementioned difference was recently used to find the direction of magnetic field in 
interstellar gas ( Heyer et al.||2008 |. 

All these techniques measure magnetic field in different media and have their limitations. For instance, the technique, 
based on the analysis of spectral lines is not applicable to hot rarefied halo gas with little emission. Studying aligned 
dust in the latter environment is also extremely challenging. At the same time, this is the domain of the new technique. 
Needless to say that obtaining the same direction of the magnetic field with different techniques increases the confidence 
of the result, while the discrepancy in the directions may testify on the existence of different regions with different 
conditions and different direction of magnetic field along the line of sight. 



6.3. Towards a cookbook for synchrotron studies 

Turbulence in our galaxy is believed to span from hundreds of parsec to hundreds of kilometers. This follows both 
from the analysis of electron density fluctuations (Armstrong et al. 1994, Chepurnov & Lazarian 2010), spectral lines 
(Larson 1981, see Lazarian 2009 and ref. therein), synchrotron fluctuations (see Cho & Lazarian 2010 and ref. therein) 
etc. With synchrotron emission one expects to measure fluctuations arising from magnetic fluctuations starting from 
the largest scale to very small scales, where the physical limitations on the small scale arise either from turbulence 
damping or Larmor radius of relativistic electrons. The practical limitations on turbulence studies will arise from both 
noise and the resolution of the data. As turbulence fluctuations decrease in amplitude corresponding to the Kolmogorov 
power spectrum, i.e. with the two dimensional measured spectrum P2 ~ K~^^/^, the small scale fluctuations are more 
subjected to the effect of noise. 

The new generation of telescopes (e.g. SKA, LOFAR) present observers with lower noise and higher resolution data 
which should enable them to reliably study synchrotron fluctuations over a wider range of scales. This will open 
prospects of detailed studies magnetic turbulence not only in Milky Way but also in supernovae remnants, radio lobes 
and external galaxies. For this wide variety of objects electron spectral index varies, but our present work shows that 
this is not a problem. 

Our present paper opens a new dimension in studies of synchrotron maps. Apart from getting the spectral index 
of magnetic turbulence, our analytical results testify that by the analysis of synchrotron anisotropics one can obtain 
the direction of magnetic fleld. Moreover, it can also deliver very unique piece of information, i.e. the relative contri- 
bution of compressible versus incompressible motions. Indeed, our study of the spatial multipoles of the synchrotron 
correlations shows that the signatures of Alfvenic and fast- wave turbulence are different. These signatures can be 
studied observationally to identify regions with different degrees of compressibility. Does higher or lower compressions 
of magnetic field correlate to higher star formation rate? What about compressibility and magnetic field amplification? 

^ We believe that this model may represent only transient state of the actual MHD turbulence and is not physically motivated in terms 
of cascades of turbulent energy. 
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These can be the questions that one wih be able to answer analyzing synchrotron maps and correlating the results 
with other astrophysical data. 

6.4. Statistics of synchrotron polarization 

While most of the present paper is devoted to studies of intensities, in Appendix E we provide the description of the 
fluctuations of other Stocks parameters and their combination. This provides an important extension of the approach 
suggested in this paper. 

We noticed that the structure of the correlations of other Stocks parameters is similar to that of intensity. On the 
basis of this we could conjecture that the dependence of the correlations involving polarization will be similar to the 
correlations of intensity, i.e. with the dependence on the cosmic ray spectral index affecting mostly the prefactor of 
obtained correlations. Numerical studies of synchrotron fluctuations with cosmic rays corresponding to different 7 in 
Ensslin et al. (2010) are suggestive of this outcome. 

In any case, in Appendix E, in analogy with the rest of the paper, we provided calculations of correlations of other 
Stocks parameters and their combinations. We observe that that combining intensity and polarization statistics one 
can remove the degeneracies, for instance, arising from the unknown 3D direction of the mean magnetic field, which 
determines the axis of anisotropy for magnetic turbulence. The bonus of combining the statistics of polarization and 
intensity could be determining of this 3D direction. This study should be continued. 

Our study opens avenues of addressing subtle issues related to magnetic field structure. For instance, Waelkens et 
al. (2009) proposed to study tension force spectrum from fourth order statistics of polarized data. They, however, 
used the isotropic model of turbulence, which does not correspond to what we now know about magnetic turbulence. 
Our approach allows such studies within a realistic model of magnetic turbulence. 

In Appendix [E] we gave, for completeness, the expressions for the sky-projected structure functions between the 
line-of-sight components of the magnetic field. While not appearing in synchrotron intensity studies, the line-of-sight 
part of the magnetic field is important in describing the Faraday rotation of the polarization direction. 

Having the statistical description of synchrotron fiuctuations and the effects of Faraday rotation on the fluctuations 
we open prospects of describing the statistics of synchrotron polarization data cubes which contain synchrotron polar- 
ization map for different frequencies (see Haverkorn 2010, Gaensler et al. 2011). Our approach opens ways to provide 
detailed quantitative description of such data cubes. 

6.5. Synchrotron foreground for CMB and high-redshift HI studies 

It is important to stress, that the description that we obtained is important for both studying galactic and extra- 
galactic magnetic fields and for removing the synchrotron fluctuations from cosmological CMB fiuctuations as well as 
from the fiuctuations high redshifted atomic hydrogen which can provide a unique insight into the processes at the 
beginning of the Universe. 

We expect our description to be important for studies of galactic synchrotr on foregrounds. In term s of weeding out 



the foreground fluctuations, a procedure of partial filtering was developed in Cho & Lazarian (2010). It was proven 
there that the galactic foreground fiuctuations, e.g. synchrotron fluctuations, arise from galactic turbulence with the 
spectral index of fiuctuations that depends on the spectral index of the underlying turbulence and the geometry of the 
emitting turbulent volume. This determinism allows one to predict the fiuctuations of foreground emission at a given 
scale, provided that the foreground fluctuations are measured at a different scale. The procedure is applicable to both 
spatial filtering of CMB signal and high redshifted HI emissi on from the early Unive rse. For instance, if synchrotron is 



separated in low resolution measurements, the procedure in Cho & Lazarian (2010) allows to remove its contribution 
for high resolution measurements. The accuracy of the removal increases if we know more about the expected power 
law of the spatial fiuctuations. In view of this the effect of 7 not being equal to the valu e 2 assumed in the estim ates 



was a concern. Our present work removes this concern and allows to use the procedure in Cho & Lazarian (2010) with 
higher confidence. 

Our obtaining in Appendix E of the tensor of synchrotron polarization should also help to weeding off the polarized 
galactic synchrotron foreground. This foreground is known to be important for the CMB polarization studies, in 
particular, for the studies of the illusive B-modes. 

In general, the more we know about the synchrotron statistics, the more reliable is its removal. For instance, if the 
filtered signal show the anisotropics expected from the synchrotron emission, this would testify that the synchrotron 
foreground was not properly removed. Additional, more sophisticated procedures may also be developed. For instance, 
the procedure of foreground removal proposed in Cho & Lazarian (2010) requires the best possible knowledge of the 
foreground spectrum at one spatial resolution in order to remove the foreground at another resolution. The a priori 
knowledge of what we expect in terms of the anisotropy of the synchrotron statistics help to determine the foreground 
statistics with higher accuracy. 

6.6. Anisotropic part of synchrotron correlation tensor 

In this paper we provided a detailed analysis of the expected properties of the anisotropics of synchrotron fluctuations. 
These anisotropics were also of interest for earlier studies, which were strongly limited in applicability because both 
of the assumption about the assumed single value of electron spectral index and of the assumed models of magnetic 
turbulence. These earlier studies are of interest as the approach presented in our paper allows to cure their deficiencies. 



For instance, a technique for obtaining the size of galactic halo was proposed in Lazarian & Chibisov ( 1991) . The 
authors proposed to use HII regions as distance indicators and calculate anisotropic part of synchrotron intensity 
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correlations for the directions towards and away from the HII region. Choosing sufficiently long radio wavelengths for 
which HII regions are opaque, one can probe magnetic fields on the scale of the distance from the observer to the HII 
region. These correlation functions were to be compared with the correlation functions obtained for both lines of sight 
avoiding HII region or with correlation functions at the wavelengths for which the HII regions are transparent. 

It is possible to see that if the angular scale of synchrotron fluctuations is 9 and the distance to the HII region is 
D, the physical size of the magnetic fluctuations is ~ DsuvO. If the angular scale of synchrotron fluctuations is 9\ 
then the ratio of the halo size to the size of the distance to the HII region will be sin 9/ sin^'. Our study of the weak 
dependence of synchrotron statistics expressions on 7 gives more confidence that similar studies using anisotropic part 
of the correlation tensor can be done in practice. 



6.7. Effects of cosmic rays 

The fact that the expression for synchrotron intensity depends on the spectral index of cosmic rays has been an 
impediment for the quantitative studies of the statistics of synchrotron variations. Our study shows that cosmic ray 
spectral variations change only the amplitude of the synchrotron variations and do n ot change the angular dependences 



of the correlation and structure functions of the synchrotron intensities (see Eq. (81 ). They do not change the direction 
of the anisotropies either. This simplifies the interpretation of the observations. 

Our study has been performed in the assumption of no magnetic field correlation with the density of relativistic 
electrons. This assumption corresponds to observation of much smoother distribution of relativistic electrons compared 
to the distribution of magnetic fields in the Milky Way. Our formalism, however, may be extended for the case of 
relativistic electrons correlated with magnetic fields. We do not expect substantial changes in our present conclusions, 
e.g. in the dependence of synchrotron intensity spectrum on 7. 



6.8. Complex geometry of observations 

In the paper we assumed that synchrotron emission is coming to the observer from a single region. Even in this 
case, the statistics of fluctuations can change from geometrical effects of observing galactic turbulence. For angles 
(f) between lines of sight less than Lturb/ Dobj where Lt^rb is th^ turbulence injection scale and Dobj is the size the 
observed region size along the line of sight, the spectral index of the angular spectrum of synchrotron fluctuations 
coincides with the spectral index of synchrotron emissivity fluctuations. For larger angles the spectral index changes. 
In an idealized case of uniform turbulence surrounding an observer, it gets universal corresponding to the correlation 
function of intensity ^ 0^^. This effect was first discussed in Shutenkov & Lazarian (1990) and attributed to the 
converging lines of sight geometry of observations. Later the change of the spectral index of synchrotron fluctuations 
was confirmed with the analysis of the observational data^" in Cho & Lazarian (2002). This paper reveals the relation 
of the spectral index of synchrotron emissivity with the spectral index of underlying magnetic turbulence which shows 
that for (j) < (j)max = Dobj/L we get the actual spectrum of the underlying magnetic turbulence. 

For studies of synchrotron emission in a volume that does not surround the observer, as this is the case for extra- 
galactic studies, it is possible to show that the direct relation of the spectral index of angular fiuctuations and the 
underlying magnetic turbulence is preserved if (p < 9max = ^obj/L, where Aobj is the extend of the emitting volume 
along the line of sight. This is the regime that one is mostly interested in. 

In terms of anisotropy studies, both regimes of > (j)max and < (pmax allow determination of the magnetic field 
direction. This is true for both galactic and extragalactic studies. 

If the observations are done looking through multiple disconnected turbulent volumes with uncorrelated direction 
of magnetic field, a decrease of the degree of anisotropy is expected. Naturally, this is not the best case for studies of 
turbulence either. Provided that for all regions along the line of sight (j) is less than (j)max, one is expected to obtain 
the same power spectral index provided that for all volumes we study, synchrotron fluctuations arise from turbulence 
in the inertial range. 



6.9. Synergy of techniques for studying turbulence spectra 

In the section above we have discussed the advantages of using synchrotron fluctuations for studies of the ISM. 
However, the technique of magnetic turbulence study that we propose in the present paper should be viewed as a 
technique complementary to other existing techniques. We believe that the most interesting questions can be answered 
if properties of turbulence in different parts of the galaxy and in different ISM phases are compared. 

In addition to the magnetic field turbulence spectrum discussed in our paper it is advantageous to study also spectra 



see 



of den sity and velocity. If the study of the former is rather straight forward from the column de n sity maps 
19991, studying velocity is far from trivial. In fact, our studies in|Esquivel & Lazarian] (|2005|; |Esquivel et al 



Stutzki 



revealed that a routinely adopted tool for turbulent velocity studies, namely, velocity centroids^^ cannot reliably reveal 
velocity fluctuations for supersonic turbulence. However, two new techniques termed Veloc ity Channel Analysis (VGA) 



and Velocity Coo rdinate Spectrum (VCS) have been introduced during the last decade ( Lazarian fc Pogosyan |2000 
2004 2006 2008 1 . These techniques have been successfully tested and applied to observational data ( Lazarian et al. 



It is worth mentioning that the change of the spectral index reported in Cho &; Lazarian (2002) on the basis of the analysis of data 
sets for high latitude synchrotron emission as well as for synchrotr on emission fr om the disk is consistent with the injection scale of 100 pc. 

One may also use so-called "modified velocity centroids" ^Laz arian fc Esquivcl 2003), but those, while increasing the range of 
applicability of centroids and increasing their accuracy in tracing velocities, still cannot be applied for the observational studies of high 
Mach number turbulence. 
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2001] iPadoan et al.|2006l [20091 iChepurnov & 



Lazariaii|2009| Chepurnov et al.|2010 see also review in Lazarian ( 2009 1 
ably retrieve data on velocity spectra from the Doppler shitted emission 



for tne references tnerein). 'i'hey allow to reliably retrieve data on velocity spectra from the DoppL 
and absorptions lines. 

We would like to stress that the overall shape of the turbulence spectrum and its spectral index can be very 
informative. The former provides information about the sources and sinks of turbulence, the latter characterizing the 
cascading of energy and, possibly, the role of plasma instabilities acting in coUisionless media. Thus advantages and 
insight arising from the ability of obtaining of magnetic field spectra and comparing its to spectra obtained in Doppler 
shifted lines are difficult to overestimate. 

Spectrum is not the only measure of turbulence. New techniques to study sonic and Alfvenic Mach numbers (i.e. 
Ms and Ma, respectively) of turbulence have been developed recently. They include so called Tsallis statistics (see 



Burlaga k F. Vifias 2004a.b, 2005. 2006., Buiiaga et al.|2006a|bl [20071 [20091 [Esquivel & Lazarian|2010, Toflicmire et al. 



2011 1, bi-spectrum (see Burkhart et al. 2009 2010), genus analysis (see Chepurnov et al.||2008 ), etc. They allow us to 
get a comprehensive picture of magnetic turbulence. Our analysis also provides a synergetic way of studying magnetic 
field compressibility. For instance, synchrotron fluctuations may map extended turbulent volume of the galactic halo. 
The relation of this turbulence taking place in tenuous plasma to the turbulence of atomic and molecular species 
which can be revealed with the VGA and VCS techniques. This will be very valuable for understanding of the complex 
interrelated dynamics of ionized and neutral media in our galaxy. 

In addition, the detailed work that we performed in this paper on describing of anisotropy of synchrotron fluctuations 
arising from anisotropic magnetic turbulence will benefit other directions of turbulence studies. Similar anisotropics 
arise, for instance in velocity fiuctuations. These anisotropics were suggested by Lazarian, Pogosyan & Esquivel (2002) 
as a means of studying the magnetic field direction and were shown to be present both in channel maps of position- 
position-velocity (PPV) data cubes of synthetic data and in the statistics obtained with velocity centroids. Later, 
velocity centroids were studied in Esquivel & Lazarian (2005, 2011) with synthetic data corresponding to numerical 
simulations with different Mach numbers, Mg and Ma, and empirical dependences of the degree of anisotropy on 
these numbers was obtained. Applying the analytical approach developed in this paper it is possible to separate 
contributions of the compressible and Alfvenic parts of the magnetic turbulence similar as we described above. This is 
going to provide new ways of detailed quantitative studies of interstellar turbulence using the wealth of spectroscopic 
surveys. 

7. SUMMARY 

Our study has revealed the following important properties of the statistics of synchrotron fluctuations: 

1. In terms of functional dependence, the spectrum of spatial synchrotron fluctuations marginally depends on the 
spectral index of relativistic electrons, while the amplitude of the fiuctuations depends on the relativistic electron 
spectral index. 

2. Magnetic turbulence makes fluctuations of synchrotron radiation anisotropic with the direction anisotropy given by 
the averaged mean magnetic field. As we found that the quadrupole anisotropy is the most prominent, this provides 
a new way of studying the direction of magnetic field using just synchrotron intensity variations. 

3. The relation between the isotropic and anisotropic parts of the synchrotron intensity fluctuations provides a way 
to study the compressibility of magnetic turbulence. 
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APPENDIX 
A. TABLE OF NOTATION 
B. AXISYMMETRIC STATISTICS OF VECTOR FIELDS 
B.l. From power spectrum to correlation function 
Let us consider axisymmetric MHD turbulence that is described by a single power E — type power spectrum 

E{k) = E{k, k • A) (^S,j - hkj^ (Bl) 

This is not the most general case of axisymmetric spectrum, since in general one can have four independent spectral 
functions. However, this ansatz arise in the regime of strong incompres sible MHD turbulence w here contribution of 



Alfvenic and slow modes is similar, as found in numerical simulations by Cho & Lazarian (20031. 

Let us consider what structure functions correspond to these models. We start from the general expression for the 
correlation function 

(i?,(xi)i7j(x2)) = J d^ke'''-"S(k) (% - fc,^,) (B2) 
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TABLE Al 
List of notations 



A unit vector of direction of the statistical symmetry axis (mean magnetic field) 

z-axis is along the line of sight 

y-axis is perpendicular to both the line of sight and the symmetry axis A 

X-axis is perpendicular to the line of sight and parallel to the sky projection of the symmetry axis A 

xi , X2 . . . 3D position vectors 

r, r 3D separation vector and its magnitude, r = (R, z) 

k, k, k 3D wave vector and its magnitude, k = (K, k^), and the unit vector k = k/fe 

fj, cosine of the angle between the 3D wave vector and the symmetry axis = k ■ A 



A 2D unit vector of the sky projection of the symmetry axis A 

Xi , X2 . . . 2D position vectors 

R, R 2D sky projection of the separation vector and its magnitude 

K, K, K 2D wave vector and its magnitude, and the unit vector K = K/K 

1/) polar angle between 2D wave vector K and the x-axis, K = (cos 1/1, simp) 

(j> polar angle between 2D separation R and the x-axis 



Isync{y^) intensity of synchrotron emission 

3y„c(R)i Daync{^) correlation/structure function of synchrotron intensity 

sy„c(R), Dsync(R) normalized correlation/structure function of synchrotron intensity 

Z)n(R) coefficient of the multipole decomposition of Dsync{K) 



Hj_(x) magnitude of the component of the magnetic field orthogonal to the line of sight. 

in'' ('")' ^h'^ i^) correlation/structure function of 7's power of Hj_{r) 

^ij-r (r), D^i (r) normalized 3D correlation/structure function of 7's power of H±{r) 



-E(k), -F(k) 3D power spectra functions. For axisymmetric turbulence -E(k) = E(k, fi) and -F(k) = F{k,ii) 

El, Fi multipole coefficient of Legendre expansion of E{k,fi) and F(k,fi) in /i 

E{ix), F{ii) angular dependence of E{k,fj,) and F{k,fi) if scale dependence is factorized 

El, Fi multipole coefficient of Legendre expansion of E{fi) and F{ii) 

cr| I , <T± variances of the magnetic field components parallel and perpendicular to the symmetry axis (mean field) 

To- ratio of (T|| and l/2a± 

-E(K), F(K.) -E^(k) , F{]<.) with kz = 0, power spectra of the fluctuations integrated along the line of sight 

En, F„, En, Fn Coefficients of the harmonic expansion of -E(K) and F(K.) or their factorized version 

cfxx , cTyy variances of the x and y magnetic field components 

e anisotropy measure of the variance of the sky components of the magnetic field, e = "^^ "yy 

Dij, Dzz 2D sky and line-of-sight components of the z-integrated correlation (structure) tensor of the magnetic field 

D+, D-, D'< D+ = Dxx + Dyy, D+ = D^x - Dyy, Dx =Dxy 



Ao unit vector of direction of the global symmetry axis in case of mean field wandering 

X average cosine of the angle between the local and the global symmetry directions 

Wi{x) weight of the isotropized spectral part 

Wl{x) weight of the local anisotropic spectral part 

G^{tp,9), G^{9) geometric functions describing the 2D structure of the turbulent mode 

Gn{S), Gn(d) multipole decomposition of the geometric functions 

rj energy injection scale 

/9 plasma constant 

Ma Alfvenic Mach number 



and rewrite it as 



(i?.(xi)i7,(x2)) = 



introduce the following differential operators 

1 d 



dr 





■ d d 








■ d d 






a 


= r- A 




1 d 



- <5., A 



dkdn]^(^^''E{k,\<i-\) 

$(r, r/i) 
r/i. Following 



(B3) 



Chandrasekhar 



d 



dr 



ld_ 
r dfx 



(1950), it is convenient to 
(B4) 



d 

a=const ' ^ 

Notice, that Dr reflects partial derivative with respect to r when a, rather than n is kept constant. is actually 
derivative with respect to a and has dimensions of 1/r. 
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In terms of these differential operators, the tensor of second partial derivatives and the Laplacian of an axisymmetric 
function are 
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dridvj 



<^{r, rn) = nrjDrr^ + SijDr^ + XiXjD^^^ + (riXj + rjXi^ Dr^^ 



and one obtains 

(7Jj(xi)Hj(x2)) = {r'^Drr^) hfj + - A) ^Sij + D^^^XiXj + rL»r^$ [fiXj + fjX^ 

that leads to the following identification of the correlation coefficients 

= r^Drr<^ , = (Z?, - A) $ , Q = D^^^ , = rl?,^* . 

B.2. Multipole decomposition 
Decomposition of the plain wave and the spectrum in spherical harmonics 

I 

e'^- = 4nY;^i'ji{kr) J2 YUmi^ir) 

I m——l 

E{k, k • A) =^ — E;(fc)y4(fc)y;„(f) 

Im 

gives after integration over angles rfOk 

$(r,r/x) =47r^z'P,(/x) fdk Mkr) Ei{k) . 
I 

The multiple expansion of the correlation tensor coefficients follows 

00 

A^ = A-kY, i'^Pnir ■ A) [-t„ + (3 + 2n)f„ - 2(1 + 2n)f„+2 + (1 + 2n)(3 + 2n)A„+4 

00 00 

+ 27r^P„(f • A)(l + 2n) ^ [4f, - + 1) - n(n + 1)) A,+2 

n=0 
00 

B^ = 47r^j"P„(f -A) 

n=0 

00 

-47r^P„(f •A)(l + 2n) ^ i'A, 



n+4,2 

+ (2n + l)A„+2 



n=0 



;=n+4,2 



Q = 27r^(l + 2n)P„(r- A) ^ - n)(Z + n + 1)A, 

i=ri+2,2 



n=0 



D^ = -2^^(l + 2n)P„(f -A) 

n=0 Z=n+1,2 



2f; - ((Z + 1)(Z + 2) - n(n + 1)) A^+a 



(B5) 
(B6) 

(B7) 
(B8) 



(B9) 
(BIO) 

(Bll) 



(B12) 



(B13) 
(B14) 
(B15) 



where T,(r) = Jk^dk ji{kr)Ei{k), f;(r) = /fc^dfc ^J±j^Eiik) and A;(r) = Jk^dk 'j^Ei{k) 

It is instructive to write out explicitly low multipole terms wrt the angle with the symmetry axis and investigate 
their behaviour as r ^ 0. Leaving only the r ^ leading terms, for the monopole 



Ao(r)~47r jk'^dk 
Bo(r)~47r j k'^dk 
Co(r)~-127r j k'^dk 



Jo(fcr-) - 3^j^ j Eo - 'i^^E, + ijj^E, 



jo{kr) - 

Mkr) 
{krf 



kr 



E'2 



and for the next order 



^2 ~ 47r / k'^dk 



j2{kr)-l 



jsjkr) 
kr 



E, + l0^:4tlE,-35'^Ee 



kr 



{kr) 



(B16) 
(B17) 
(B18) 

(B19) 
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B2 - -47r J k-'dk 
C2 ~ 1407r y" k'^dk 
Di^Un [ k^dk 



[krY 



^2 + 5^i?4 



£^4 



j4,{kr) 



(B20) 
(B21) 
(B22) 



At zero separation r ^ 0, Ao(0) and L>o(0) while Bo{0) = 4tt Jk'^dk [^Eq + ^£2] and Co(0) = -f^ Jk'^dkE2 
describe the variance and same-point correlation of the components of the vector field 

(iJ,(xi)ff,(xi)) = B(0)(5y + C{0)\,\j (B23) 

Hence, for axisymmetric field not only the variances of the components differ, but there is zero separation correlation 
between different components if no frame axis is aligned with the symmetry axis. 

In case of structure functions one computes -B(O) — B{r), . . ., which leads to regularization of the monopole terms 
Bq and Co- 

B.3. Derivation of the multipole decomposition 

Derivations of the expressions Eqs. ( B12||B15 ) involves evaluating derivatives of the ji{kr)Pi[^) basis functions and 
rearranging the series in multipole "1" . We start with the following relations for the the angular derivatives 



d 



ii^^mi)=imi)+ J2 (2" + i)^n(M) 



n=(0,l),2 
l~2 



Pi 



(/i) = ?2Pi(M) + J E [(^ + 2)a-l)-n(n + l)](2n + l)/^„(M) 



n=(0,l),2 
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(B24) 
(B25) 
(B26) 



" n=(0,l),2 

where the notation X]ti=(o 1) 2 Hi^ans that the sum over n starts from the first value in the parenthesis (zero) when I is 
even and from the second one (unity) when I is odd, and proceeds with the step of two. The sums are zero for / = 0, 1. 
For the radial derivative we shall use 



1 d I 



(B27) 



Now we can write out the actions of several operators that we need. The following expressions show some intermediate 
steps, and for the final conclusion point to the structure that the operator creates in the full multipole expansion after 
the integration over k has been carried out 

i/(fcr) a _ ,2j/+i(fcr) ,2ji{kr) 



1 d 



Dr [ji{kr)Pi{^i)] = P,(^)- j;(fcr 

r or r 



2 Ma^^'(M) 



kr 



^2^2 



J2 {2n + l)P^{^,) 



Tl=(0,l):2 



r^lD^Dr[Jl{kr)Pl{^i)] = 



-tl{r)Pli^l) - Ai{r) E (2n + l)P„(/i) 

n=(0,l),2 

d 



^ig^JDr[jl{kr)Pl{^l)] 



(B28) 
(B29) 



1-2 



D^^[ji{kr)Pi{y)] 



lPi{fi)+ Y.{2n + l)P„{fi) 
"=(0,1), 2 

— i~ J \ III — ^ 



1-2 



\Ki{r) E [(^ + 2) -n(n + l)](2n + l)P„(A.) 

n=(0,l),2 

Pl{^l)^\^l{r) E {l~n){l + n + l){2n + l)P,{p) , 



(B30) 



n=(0,l),2 



^{3i{kr)PM) = -eji{kr)Pi{^i) ^ -T;(r)P,(^) 



(B31) 

Multipole expansion of the correlation coefficients now follows from combination of the above results. In order to 
obtain the final expressions, we make two important steps. Firstly we change the order of summation X]/^o X]i=(o 1 



(0,1):2 



Sn=o S/=n+2 2 ^^"^ ts^G the label for the exterior sum to be n everywhere. Secondly we collect together the lowest 
order terms of the same power of divergence in k integrals. 
Let us proceed in the order of difficulty 
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Cj is obtained directly from Eq. (B30). 



B.3.1. Q 



B.3.2. Bf 



is a combination of Eq. (B28l and Eq. (B31). Let us illustrate the reversal of the summation order on this 
example. Using n for the exterior index the direct combination of two equations and the change of the summation 
order in the double sum gives 

oo oo oo 

=47r^z"P„(^) (t„-f„) -47r^(27i+l)P„(/i) ^ i'Ai (B32) 

n=0 n=0 ;=n+2,2 

It is important that Aq contribution is not present since the integral in A„ c>c jn{kr) / [kr)"^ ~ (fer)""^ is potentially 
the most divergent one. T„ ~ jn{kr) ^ (kr)"^ and F ~ jn+i/{kr) ~ (/cr)" has the same level of potential divergence as 
A„+2- These considerations are important for the steep spectra that play role in turbulence studies, where correlation 
functions are not formally determined at long wavelengths, and transition to regularized structure function must be 
performed. Such transition amounts to regularization of the integrals, mainly in A„. Thus, we group the I = n + 2 
term from the double sum together with the first term to get 



B^^4nJ2i"Pn{t^) (T„-f„ + (2n+l)A 



n+2 



n=0 i=n+4,2 



(B33) 



B.3.3. 



Using Eq. (B6| w e represen t the coefRcient as = (A — 3Dr — Df^^ — 2r^DrDfj) $, thus reducing it to the 
combination ot Eqs. ( B28|B3r| ). Let us collect separately the terms in the single sum and the terms in the double sum 



(3 + 2n)f„ 



+ 



4^^(2n + l)P„(Ai) [2f, + i((Z-l)(Z-2)-n(n + l))A, 



(B34) 

n=0 /=n+2,2 

We again separate the first term in the double sum, noting that there is an additional cancellation of A„_|_2 contribution 



A^^A^Y. + (2n + 3)f„ - 2(2n + l)f„ 



n=0 



+2 



4^^(2n + l)P„(Ai) [2f, + i((/-l)(Z-2)-n(n + l))A, 



(B35) 

n=0 /=n+4,2 

and then shift ^ — Z + 2 in the A part of the double sum to select A„+4 which has the same scaling at fcr — ^ as the 
terms in the single sum 



A5=47r^i"P„(/i) 



oo 



-T„ + (2n + 3)r„ - 2(2n + l)r„+2 + {2n + l)(2n + 3)A„+4 
2fi + \{{l + l)l-n{n+l))Ki+2 



+ 



(B36) 



+ 47r^(2ji + l)P„(/i) Y ^' 

n=0 ;=n+4,2 

which gives our original expression. 

B.3.4. 

The last coefficient, Z?j is different from the ones that have just been considered in that its multipoles are linked to 
multipoles of the power spectrum of the opposite parity. While even multipoles of the only expected to be present, 
give rise to even multipoles of Aj, and C^, they create odd multipoles of the D^. As a consequence, cannot be 
assembled using just even operators that we have studied. The additional odd derivative that we need is just 

Y ^"(^) (B37) 

^ n=(l,0),2 

(here summation starts with unity if 1 is even, and zero if 1 is odd, I — term is absent) which allows us to write 

1-2 



rD^Dr {ji{kr)Pi{p)) = ^ [I?,ji(fcr)P,(/i)] 



-f,(r)|-P,(/i)-A,(r) Y (2n + l)|^P„(/i) 

^ n=(0,l),2 " 



i-1 ^ l-l 

tiir) ^ (2n + l)P„(Ai)--A,(r) ^ (/ + n)(/-n-l)(2n+l)P„(^) (B38) 



ra=(l,0),2 



n=(l,0),2 
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Reverting the summation order again in the fuh niuhipole expansion, X];^o o) 2 ~ "^nLo Sz^n+i 2' ^^'^ shifting 



the summation index in A/ part by two, Z — > ^ + 2, we obtain the Eq. (B15) 



C. CORRELATION OF THE SYNCHROTRON INTENSITY IN AXISYMMETRIC MAGNETIC FIELD 

FOR 7 2 

For 7 = 2 we have the following relations 

{Hl) = {Hl^Hl)^a.,^+Oyy (CI) 
{Hi) = {{Hl + Hlf) = 3aL + ioly + 2a^^ayy + 4a^^ (C2) 



2\2\ 

y' 

{Hi (xi) Hi (X2)) ^{{Hl + Hl)^{Hl + H^h) 



(C3) 



: {Hi (xi) Hi (X2)) + {Hi (xi) Hi (X2)) + (i^^ (xi) Hi (x^)) + (i/.^ (xi) Hi (x^)) 



[iji(xi)-i?i(x2) 



-2{H, (xi) i/, (X2))' + 2 (i/, (xi) (X2))' + 2 (i/^, (xi) iJ, (xa))' + 2 (i/j, (xi) Hy (xa))' 
= 2(i7l)-2(i7i(xi)i/i(x2)) = 



r2 



j2 . 



r2 

'xjf 



- 4 (i/, (xi) i/, (X2))' - 4 (if, (xi) Hy (X2)>' - 4 (i/^ (xi) H, (xa))' - 4 {Hy (xi) ff, (xa))' 



8<TxyDxy - 2D^y 



(C4) 

where we have used the fact that the correlation tensor is symmetric and have denoted the same point 2D covariance 
matrix by ciij = (HiHj). The final expression can be written in the manifestly 2D-invariant tensor form 



{Hi)-{Hir^2j2a.,a 
2 

(i/i(xi)i?i(x2))-(i/i)^-2^e.^- 



2 



{[Hi (xi) - Hi (X2)] ) - ^ (4a,,Z?,, - A,^j^) 
The correlation function of the synchrotron intensity signal / (X) cx J dzHl (X, z) is then given by 

(/ (Xi) / (X2)) - (/ (Xi))' '^j'^^^j '^^^ (Xi, ^1) Hi (X2, Z2)) - {Hi (Xi, zi)) {Hi (Xi, ^2))] 

= 2 y"dz+ J dz {X, z) {X, z) . 
Correspondingly, the structure function is 

I dziHl (Xi, zi) - I dz2Hl (X2, Z2) 



(C5) 
(C6) 
(C7) 



([/(Xi)-/(X2)]0« 



(C8) 



(C9) 



dz+ J dz {4a,, [i?,,: (X, z) - D,, (0, z)] - [D,, (X, z) D,, {X, z) - D,, (0, z) D,, (0, z)]} . 

D. GEOMETRICAL FUNCTIONS 

The synchrotron correlation function in approximation of the current paper is expressed via the trace of the 2D 
projection of the spectral tensor that describes the given type of waves. We have encountered the multipole expansions 
of the following angular functions that represent the 2D angular dependence of the traces of the Fast, Alfven and Slow 
projected spectra 



27r Jq 1 — cos-^ip sm 6 



1 

2^ 



dipe 



1 — cos^tjj sin 



(Dl) 
(D2) 



The two angular functions are related by Gp{d) + sin 9Gp (9) — Spo- In Figure Dl we plot the behaviour of the 
multipoles up to p = 6 as the functions of the symmetry axis orientation angle 9. 
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Fig. D1. — Low multipole behaviour for multipole decomposition of the traces of the spectral tensors for the Fast and Slow (left, general 
F-type tensor) and Alfven (right) waves. 

E. CORRELATIONS OF THE POLARIZATION COMPONENTS OF THE SYNCHROTRON 

RADIATION 

E.L Correlation of Stokes parameters 

Synchrotron intensity is just one of the characteristics of synchrotron emission. Studies of the statistics of synchrotron 
polarization can be very promising (see Waelkens et al 2009). Below we discuss how the approach in the main part of 
the paper be generalized to study other Stocks parameters and their combinations. In what follows we use the realistic 
model of MHD turbulence that we advocate in the paper. For the sake of simplicity we assume that 7 = 2. A study 
of 7 dependence will be done in the subsequent papers. 

Complete description of the radiation incoming from the direction X on the sky is encoded in the polarization matrix 



P(X) = 



u 



Q U 
iV I 



iV 



= /(X) 



1 
1 



Q(X) UiX)+iV{X)\ 
C/(X)-iT/(X) -Q(X) J 



(El) 



where Q,U,V are well-known Stokes parameters. 

Here we shall limit our consideration to the linear polarization, described by Q and U . We are interested in 
the correlation between two matrixes (P(X)P(X')). Such correlation is convenient to separate into parts based on 
transformation rules under rotation of the coordinate frame. This gives rise to the 



scalar : (/(X)/(X')) 



"vector" 



(/(X)O(X')) 

(/(X)C/(X')) 



"tensor" : 



( (Q(x)g(x')) (g(x)c/(x')) 
I (c/(x)g(x')) (f/(x)[/(x')) 



(E2) 



parts. Here we put "vector" and "tensor" in the quotation marks, since they rotate with twice the angle $ of the 
rotation of the coordinate system 

' Q' \ _ ( cos 2$ sin2$ \ / Q 
U' ) ~ \ -sin 2$ cos 2$ ) \ U 



(E3) 



The "vector" terms describe cross-correlation between the intensity and linear polarization, while the "tensor" part 

describes correlation between polarizations themselves. 

For synchrotron radiation in 7 = 2 case, following Waelkens et al. ( |2009[ ), 



2 

XX 



U{X.)(xp rie / 



(E4) 
(E5) 



(E6) 



where p — \/ Q"^ + U"^ / 1 is the degree of linear polarization and rig is the density of the electrons, both assumed 
constant over the line-of-sight range where the signal forms. These expressions are in the laboratory frame x^y in 
which Q and U parameters are defined. 



Evaluation of the structure functions that involve polarization are similar to the computations that led to Eq. ( 48 1 
for intensity. Using for brevity notation £)/[/(R) — ([/(Xi) — /(X2)] [J7(Xi) — C/(X2)]) and so on, the expressions for 
individual correlation components are 



DjQ cx 4p / dz X 



(E7) 



CTxxDxxi^, Z) - O-yyDyyCR, ^) ~ ^ (L'x:e(R-, ^)^ + DyyiH, z)^ 



1 







Diu (xAp j dz X 

^xy iPxx{^-, z) + £'to(R, z)) 4- {px 



(E8) 



Oyy) £'j;y(R, z) - -Dxy(R, z) {D.j:x{K, z) + ^^^(R, z)) 



R^ 



30 

Dqq oc I dzx 



(^XxDxxi'i^, Z) + (JyyDyyill, Z) - 2(J^yD^y{Ii, ^ (^^^^(R, + Dyy{^, - 2 D ^ y {R , 



Duu 4p dz X 



(E9) 
R-> 
(ElO) 







i^Qtr oc dz X 



CTxy {D^J:,,{R, z) - Dyy{R, Z)) + {a^^ - (Tyy) D^y{R, z) " -D^y{R, z) {D^.^{R, z) - Dyy{R, z)) 



(Ell) 
R^ 



where [R — ?> 0] signifies the same integrand as in the main expression but evaluated at zero separation R between the 
hne of sights. We remark that U or Q correlations with intensity / are proportional to the first power of the degree 
of polarization, while polarization cross-correlations are quadratic. For weakly polarized signal it may be easier to 
measure the former. 

The trace of the polarization correlation "tensor" provides the frame invariant measure 



Dqq + Duu « 4p dz 



1 2 

{(^XX + (Jyy){Dxx{R, Z) + Dyy{R, Z)) - - {Dxx{R, Z) + Dyy{R, z)) 



R^O 



(E12) 



In Eqs ( E7|E11 ) the 2D components of the magnetic field correlation (structure) tensor Dij refer to the laboratory 
frame, in which tfie Stokes Q and U are measured. Let us rewrite the same expressions using Dij in the frame aligned 

with the symmetry axis A, where theoretical calculations for axisymmetric turbulence are more straightforward. From 
here onwards, D~^ = D^x + Dyy, D~ = D^^ — Dyy and D^ = D^y refer to components in such a symmetry frame, while 

$ is the angle between A and the laboratory x-axis. Limiting ourselves to the linearized approximation, applicable, as 
in Eq. (48), at small scales, we find 



DiQ(x2p (cT^^+ayy) J dz [{D~ (R, z) + eD+ {R, z)) cos 2$ - 2L>'' (R, z) sin 2$] - [R ^ 0] 
Diu(x2p {cTxx+fTyy) j dz [(Z?" (R, z) + 61)+ (R, z)) sin 2$ + 2D''(R,z)cos2$] - [R ^ 0] 
DQQ(x2p^{(j^^ + (jyy) j dz [D+{R,z) + eD^{R,z)cos^^ - 2ei:''' (R, z) sin4$] - [R ^ 0] 
Duuo^2p^{a^^ + (Tyy) j dz [i:>+(R,z) - ei:>"(R, z) cos4$ + 2ei:>'' (R, z) sin4$] - [R ^ 0] 
Dquo^2p^ e{(J:^^+(jyy) j dz [/^"(R, z) sin4$ + 2L»''(R, z) cos4$] - [R ^ 0] 



The frame invariant combinations are 



Dqq + Duu « 4p2 (a^^ + ayy) j dzD+{R, z) - [R ^ 0] , 



and, for completeness, the intensity structure function from Eq (48) 



Dii cx 2 (d,, + ayy) J dz [{D+{R, z) + eD~{R, z))] - [R ^ 0] 



(E13) 
(E14) 
(E15) 
(E16) 
(E17) 

(E18) 
(E19) 



The magnetic field variance anisotropy parameter, e, is defined in Eq. (50). 

E.2. Angular properties of the polarization correlations 

Let us turn to the angular dependence of the polarization functions. Our study of the intensity of the synchrotron in 
the main text can be reformulated as the following rules of correspondence between the different parts of the magnetic 
correlation tensor and the multipoles of the synchrotron correlations they contribute to 



/oo 

dzD+(R,z) ^ DnOC EsGn-siO) 

oo ^ 

dzZ?-(R,z) ^„CX ^ ~-{Es+2+Es-2)Gn-s{e) 



(E20) 
(E21) 
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where Es are coefRcients of multipole expansion of the projected power spectrum E{K, ^p) wrt the sky orientation 
angle of the wave vector tp and Gp{0) are multipoles of one of the geometric functions G{9,'il)) (see Appendix IdI) that 
describes the specific turbulent mode. This rule reproduce s E qs. (61) and (65) when applied, correspondingly^ to E 
and i^-type modes of the general axisymmetric tensor Eq. 



(511 



To study polarization one also needs the additional component 



dzD^'in, z) 



oo 

y - 

s — — oo 



{Es+2 — Es-2) Gn~s{0) ■ 



(E22) 



While Z?^ and D combinations give rise to real £)„ coefRcients and, thus, symmetric ~ cos n(j) dependence of the 
synchrotron correlations, contributes the imaginary part to Z)„ and antisymmetric ^ smn(f> signal correlation 
behaviour. Antisymmetric contribution is not present if the mode tensor is isotropic, when Gn-p = Snp- Presence of 
the antisymmetric, sinn0, behaviour in the correlation patterns of polarization is a direct indication of the i<"-term 
contribution to the spectral tensor structure. It is expected at some level for Alfven and Fast modes. 
The resulting multipole decomposition of the polarization structure functions is 

-I iE,+2 + E,_2) + eE, 



p,Q)„(x2p^C„(m)i?i+™ J2 



cos 2$ + - {E,+2 - E,_2) sin 2$ 



g: 



{A,F,S) 



(E23) 



l+m 



E 



s+2 



eE, sin 2$ 



[Es 



cos 2$ 



g: 



{A,F,S) 



(E24) 



(DQQ)„(x2p2 AG„(m)i?i+™ 



E, 



2 {Es+2 + Es 



I cos 4$ 



^ {Es+2 



Es-2) sin 4$ 



G 



{A,F,S) 



(E25) 



(-Dc/c/)„cx2p2 ACn{m)R 



{Dqu)^2p^ AG„(m)i?H 



where G„_'s ' {9) stands for one of the particular geometrical mode functions, Alfven, Fast or Slow, or their combi- 
nations. It depends, parametrically on the angle 9 between the mean magnetic field direction and the line of sight. 

Polarization provides new avenues to link properties of the magnetized turbulence to observables. The most direct 
measure of the turbulence spectral anisotropy is provided by the trace of the synchrotron polarization matrix 



l+m 



E 



Es 



{E. 



s+2 ■ 



Es^2) COS 4$ - - {Es+2 - Es-2) sin 4$ 



G 



{A,F,S) 



00 

E 

S — — OC 



- - {Es+2 + Es-2) sin 4$ - - {Es+2 ~ Es-2) cos 4$ 



G 



{A.F,S) 



(E26) 



(E27) 



(X 4p^ A(^A,F,S)Gr, 



{m)R 



l+m 



E 



EsG. 



{A,F,S) 



{0) 



(E28) 



It may be possible to determine the direction of the axis of symmetry (mean magnetic field) by measuring separately 
symmetric (cos) and antisymmetric (sin) patterns of the correlations. For example tan 2$ = Re(£'/[/)„/Re(_D/Q)„ can 
be compared with the direction obtained from elongation of iso-correlation contours of intensity. And when symmetry 
direction is determined, e parameter can be estimated from Im(£'Q[/)„/Im(Z)/(7)„ — pe evaluated in the symmetry 
frame $ = 0, if there is a non- vanishing imaginary contribution to correlation multipoles from Dxy 

W e may con clude that angular and scaling dependencies of the correlations of synchrotron polarization and intensity, 
Eqs ( |E13|E19 ), contain wealth of information about orientation of the symmetry axis, e parameter and full projected 
correlation tensor of the magnetic field. Redundancies build into different signals allow for multiple cross-checks. 
Detail study of the possibilities that synchrotron polarization opens for studies of anisotropic magnetized turbulence 
is a subject of the future work. 

E.3. Correlations of the line-of-sight component of the magnetic field 

Faraday rotation measures provide a valuable source information about magnetic turbulence. However, the models of 
turbulence adopted in earlier studies were models of isotropic turbulence. Below we calculate the statistics of Faraday 
rotation measures which is based on a realistic MHD turbulence model. The measures we deal with do not depend on 
the cosmic ray spectral index. 

Study of intensity and linear polarization of the synchrotron led us to the discussion of the orthogonal to the line- 
of-sight components of the magnetic field. Let us in conclusion collect the results for the projected correlation of 
the line-of-sight component /^^^(R) — J dz{Hz{:x.i)Hz{x2) ■ This quantity describes the correlation of Faraday phase 
((/)(Xi)0(X2)) and is important in polarization studies. We present the results for the line-of-sight component here to 
complete ou r st atistical formalism. 

From Eq. (51 1 one immediately obtains the sky-projected structure function between the line-of-sight (z) components 



of the magnetic field as 



I?,,(R) - / dz{H,H, 



1 



(2^)2 



E{K) + F{K) 



cos^9 



1 ~ cos^V sin 



(E29) 
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Spectral representation of different turbulence modes ca n be expressed again using combinations of the geometrical 
functions introduced in Appendix [d| Namely, using Eqs. 



Alfven 

Slow (high /?) 
Alfven + Slow (strong) 
Fast [low /3) 
Fast {high P) 



(I69|[77 

F(K) - -E{K) 

E{K) = 0, 

F{K.) = 
E{K) = 0, F{K) = F{K) 
E{K) = 0, ^^(K) = F{K) (1 



79 73) 



cos^Vsin^6') Dz 



cx ^;(K) sin^eiG^ 
cx F{K)G^ 
o^E{K) 
cx F{K)G^ 
(X F{K) cos^e 



(E30) 



Multiple decomposition of the angular dependence of the D^z (R) is described by convolution of the multiple decom- 
positions of the spectral and geometrical functions in Eq. (E30). Variation of the direction of the mean magnetic field 
along the line of sight leads to isotropization of the 
intensity. 



m a similar way that we have studied for the synchrotron 
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